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SUMMARY 


This  report  presents  the  results  of  a  research  and  development  effort  performed 
under  Air  Force  Contract  F33615-77-C-3121,  entitled:  "Improved  Methods  for  Pre¬ 
dicting  Spectrum  Loading  Effects. "  The  objective  of  this  program  was  to  update  the 
crack-growth  prediction  technology  required  for  implementation  of  the  damage- 
tolerance  control  procedures  throughout  the  life  cycle  of  any  weapon  system.  Pri¬ 
mary  efforts  devoted  to  development  of  a  preliminary  design  level  damage-tolerance 
analysis  method  were  first  to  formulate  a  crack-growth  life-prediction  method  suitable 
for  use  in  the  preliminary  design  stage  of  configuration  development.  A  crack-growth 
prediction  module,  PREGRO,  which  utilizes  a  spectrum  characterization  method,  was 
then  incorporated  into  the  Automated  Pre-design  of  Aircraft  Structure  (APAS  III),  a 
multistation  structural  synthesis  procedure  of  the  Structure  Technology  Evaluation 
Program  (STEP).  The  figure  on  page  vi  is  representative  of  the  overall  STEP 
system.  The  General  Interactive  Executive  Management  System  (GEMS)  is  the 
executive.  Part  of  GEMS  is  an  application  executive  (APEX).  The  Application 
Data  Manager  (ADM)  provides  the  data  base  management  and  pre/post  processing 
of  data  for  application  programs.  APPL  shown  in  the  figure  represents  an 
application  program  (batch)  with  an  APPLIN  input  file  and  APPLOUT  output  file. 
APPLGR  represents  an  interactive  graphics  application  program.  The  figure  on 
page  vii  shows  the  relationship  of  the  various  application  program  modules  to 
the  data  files  preprocessed  by  the  Application  Data  Manager  (ADM)  and  output 
data  files  post -processed  by  the  ADM  for  inclusion  in  the  data  base. 

PREGRO  performs  a  cycle -by -cycle  crack- growth  analysis  of  a  uniblock  flight 
spectrum  which  accounts  for  the  load  interaction  effects  to  obtain  the  crack-_ 
growth  rate  per  flight  (da/dF)j  and  a  measure  of  the  stress  intensity  factor  Kj 
for  j  values  of  initial  crack  sizes.  It_then  characterizes  an  equivalent  growth- 
per- flight  equation  (i.e. ,  da/dF  =  ck*)  to  obtain  the  growth  rate  parameters  c, 

\  ,  and  ( Aa2)  1/2  through  a  least -square -fit  routine.  The  final  crack-growth 
life  is  calculated  through  the  linear  approximation  technique. 

Supplementary  effort  was  directed  toward  extending  the  utility  of  APAS  for  conducting 
rational  preliminary  design  trade  studies  on  different  aircraft  categories.  These  tasks 
consist  of  restoration  of  the  full  range  of  stress  intensity  factor  correction  function  for 
riveted  stiffened  panels  containing  large  cracks,  an  extension  of  the  crack-growth 
analysis  process  to  include  unstiffened  plate  construction  concepts,  an  extension  of  the 
load  spectrum  library  to  include  a  typical  spectrum  for  a  lightweight  fighter  aircraft, 
and  incorporation  of  an  alternate  means  for  specifying  load  spectrum. 
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This  report  was  prepared  by  the  Rockwell  International  Corporation,  North 
American  Aircraft  Operation,  Los  Angeles,  California,  under  Air  Force  Contract 
F33615-77-C-3121,  Project  2401  "Structural  Mechanics,"  Task  240101,  "Structural 
Integrity  for  Military  Aerospace  Vehicles,"  Work  Unit  24010120,"  Improved 
Methods  for  Predicting  Spectrum  Loading  Effects . "  This  work  was  administered 
under  the  direction  of  the  Air  Force  Wright  Aeronautical  Laboratory,  Flight 
Dynamics  Laboratory,  Structure  and  Dynamics  Division,  Wright -Patterson  Air 
Force  Base,  Ohio.  Mr.  Robert  M.  Engle  was  the  project  engineer. 

This  work  was  accomplished  by  personnel  from  the  Fatigue  and  Fracture  Mechanics 
Group,  Dynamics  Technology,  Structural  Systems,  supervised  by  George  E.  Fitch,  Jr. , 
supervisor,  Joseph  S.  Rosenthal,  manager,  and  Dr.  Leslie  M.  Lackman,  director. 
James  B.  Chang  was  the  program  manager  and  principal  investigator.  The  major  task 
performed  was  to  implement  a  crack-growth  analysis  procedure  which  accounts  for  the 
load  interaction  effect  of  the  flight  spectrum  to  crack  growth  into  the  APAS  in,  a 
structural  synthesis  procedure. 
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SECTION  I 


INTRODUCTION 


Automated  Pre- Design  of  'Aircraft  Structure  (APAS)  is  a  multistation  structural 
synthesis  procedure  designed  to  perform  structural  sizing  and  detail  design  box 
beam  structure.  Development  of  this  program  was  started  at  General  Dynamics 
Convair  Division  in  1972.  The  original  version  of  APAS  was  limited  to  static- 
strength  considerations.  The  procedure  was  further  developed  to  provide  the 
addition  of  fatigue  and  fracture  mechanics  design  criteria  and  multimaterial 
capability  including  advanced  composite  materials.  This  version  of  APAS  was 
identified  as  APAS  III ,  which  represents  a  unique  capability  for  evaluating  the 
weight  impact  of  various  design  criteria  that  may  be  considered  for  the  primary 
structure  (fuselage  shells  and  wing  boxes)  of  transport  aircraft.  Criteria  that 
may  be  evaluated  by  APAS  III  include  material  selection ,  static  loads ,  geometry , 
structural  configuration,  minimum  gage  limits,  fatigue  life,  crack- growth  life, 
and  residual  strength.  Reference  1  described  the  APAS  III  program  in  detail. 

The  crack-growth  module  in  APAS  III  lacks  the  capability  of  accounting  for  load 
interaction  effects  such  as  tensile  overload  retardation  and  compressive  load 
acceleration.  Therefore,  to  utilize  APAS  HI  in  the  preliminary  design  trade-off 
studies,  a  false  indication  on  the  calculated  weight  panelty  might  occur  when  a 
candidate  structural  component  is  sized  based  on  the  crack-growth  analysis  result. 
Under  Air  Force  Contract  F33615-77-C-3121,  "Improved  Methods  for  Predicting 
Spectrum  Loading  Effects,  ”  a  research  effort  has  been  devoted  to  develop  a 
crack-growth  module,  PREGRO,  which  realistically  accounts  for  the  load  interaction 
effects.  This  crack-growth  module  was  then  incorporated  into  a  modified 
APAS  III.  The  modified  version  of  APAS  III  is  identified  as  APAS  IV  in  this 
report . 

Supplementary  effort  was  also  directed  toward  several  tasks  which  are  aimed  at 
extending  the  utility  of  APAS  for  conducting  rational  preliminary  design  trade 
studies  on  different  aircraft  categories.  These  tasks  consisted  of  the  restoration 
of  the  full  range  of  stress  intensity  factor  correction  functions  developed  by  Poe 
for  riveted  stiffened  panels  containing  cracks  into  APAS  IV ,  extension  of  the 
APAS  IV  crack- growth  analysis  process  to  include  unstiffened  plate  construction 
concepts ,  extension  of  the  load  spectrum  library  of  APAS  IV  to  include  a  typical 
spectrum  for  a  light-weight  air-to-air  fighter,  and  incorporation  of  an  alternate 
means  for  specifying  load  spectra.  Reference  2  documents  the  development 
efforts  of  the  abovementioned  tasks.  This  report  provides  the  guidance  to  the 
user  for  executing  APAS  IV. 


SECTION  II 


COMPUTER  PROGRAM  CAPABILITIES 

2 . 1  GENERAL  DESCRIPTION 

The  technical  approach  used  in  APAS  IV  is  applicable  to  any  closed  section  beam-like 
structure,  and  it  is  typical  of  the  procedure  used  in  the  early  design  phase  of  aircraft 
structure.  The  overall  approach  makes  use  of  a  point  design/analysis/redesign  pro¬ 
cess  that  is  iterated  until  an  acceptable  design  is  produced.  Figure  2-1  presents  the 
functional  flow  chart  for  the  approach  used  in  the  APAS  FV  computer  program.  This 
flow  chart  outlines  the  major  analysis  and  design  loops  of  the  program.  A  detailed 
discussion  of  the  theoretical  basis  for  the  analysis  is  presented  in  Section  4.1.3, 
LIMITATIONS  AND  RESTRICTIONS. 

2.2  FUNCTIONS  PERFORMED 

The  APAS  routines  perform  structural  sizing  and  detail  design  of  box-beam  structure. 
Inputs  to  the  program  include  component  geometry,  structural  elements,  flight  pro¬ 
file/load  spectrum,  and  external  loads.  The  structural  analysis  may  select  from  the 
following:  criteria,  materials,  geometry,  static  strength,  stability,  fatigue,  fracture, 
and  residual  strength. 
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SECTION  in 

FUNCTION  DESCRIPTION 


3. 1  TITLE  OF  FUNCTION 

This  section  discusses  APAS,  the  automated  predesign  of  aircraft  structure  program. 

3.2  DESCRIPTION  OF  FUNCTION 

The  APAS  IV  program  is  highly  modular  to  facilitate  modifications  and  additions. 

This  modular  approach  has  resulted  in  a  large  number  of  subroutines  which  are  briefly 
described  below.  A  complete  layout  of  functional  relationships  of  these  subroutines 
is  presented  in  Figure  3-1. 

A  skeletal  framework  was  initially  laid  out  for  APAS  IV  with  dummy  subroutines  to 
provide  for  future  expansion  of  the  program.  Some  of  the  dummy  subroutines  still 
exist,  and  the  functions  meant  to  be  performed  by  these  routines  are  also  described 
below.  The  dummy  routines  mainly  concern  themselves  with  additional  input  and 
output  options  and  do  not  affect  the  basic  problem  solving  functions  of  the  program. 

APAS 

This  is  the  overall  driver  program  that  calls  the  input,  output  and  processing 
routines. 

BOXLDS 

This  subroutine  performs  a  box  beam  internal  load  solution  at  a  cross  section 

and  computes  the  complex  bending  stresses  and  shear  flows  for  unit  load  components. 

BUCKLE 

This  subroutine  performs  a  panel  buckling  analysis  for  general  instability  of  simply 
supported  curved  orthotropic  panels  in  bi-axial  loading  with  shear. 

CHKBT 

This  subroutine  checks  the  input  BEE  and  TEE  variables  for  illegal  values. 

CLAMDA 

This  subroutine  performs  a  least  square  fit  procedure  to  calculate  the  power 
exponent ,  ’P  and  the  growth  rate  constant ,  c ,  for  an  equivalent  crack  growth 
per  flight  rate  equation. 

CONVRG 

This  subroutine  determines  when  to  stop  the  redesign  iteration  process  based  on 
the  input  convergence  criteria  parameters. 

CORREC 

This  subroutine  calculates  the  stress  intensity  factor  width  correction  factor  and 
combines  it  with  the  other  stress  intensity  correction  factors. 
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CPLIM 


This  subroutine  is  a  recovery  procedure  which  will  be  executed  in  the  event  of  a 
CP  time  limit.  It  causes  the  values  of  output  variables  at  the  time  of  CP  time 
limit  to  be  output.  This  subroutine  will  be  nonfunctional  if  the  subroutine  RECOVR 
is  not  available. 

C  RIFLE 

This  subroutine  calculates  the  crippling  stress  of  a  panel  stiffener  using  non- 
dimensional  crippling  curves. 

CRIPRP 

This  subroutine  calculates  the  section  properties  of  the  stiffened  panel  configura¬ 
tions,  including  effective  skin  if  flagged. 

CRITIC 

This  subroutine  solves  for  the  critical  crack  length  for  an  infinitely  wide  plate. 
CRKDAT 

This  subroutine  fills  the  stress  intensity  factor  array  and  the  stiffener  load  con¬ 
centration  factor  array  for  the  case  of  riveted  stiffeners. 

CRKSIZ 

This  subroutine  fills  the  stress  intensity  factor  correction  factor  array  for  the 
case  of  integral  stiffeners. 

CUBE 

This  subroutine  fits  a  cubic  polynomial  to  four  points  and  finds  the  minimum 
value.  The  minimum  found  is  compared  with  the  minimum  found  on  the  previous 
call.  If  they  agree  within  tolerance  then  the  convergence  flag  is  set.  This  routine 
is  used  by  the  one  dimension  search  subroutine  ONED. 

EFFSKN 

This  subroutine  computes  the  effective  width  of  skin  acting  with  a  stiffener.  The 
skin  and  stiffener  materials  may  be  different. 

EVA 

This  subroutine  evaluates  the  overall  acceptability  of  a  structural  element  based 
on  manufacturing  constraints  and  stress  analysis. 
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EVAG 


This  subroutine  is  called  by  REDSON.  It  determines  the  flaw  growth  margin  of 
safety  for  each  element  in  a  symmetry  group. 

EVAL 

This  subroutine  returns  the  margin  of  safety  for  fatigue. 

EVAR 

This  subroutine  is  called  by  REDSON.  It  determines  the  residual  strength  margin 
of  safety  for  each  element  in  a  symmetry  group. 

FATTB1,  FATTB2,  FATTB3 

These  BLOCK  DATA  subroutines  store  fatigue  data  in  the  form  of  constant-life 
diagrams. 

FATTIN 

This  program  reads  a  set  of  S-N  data  from  input  cards  or  loads  a  set  from 
FATTAB. 

FLTGRO 

This  subroutine  performs  the  integration  of  the  crack  growth  per  flight  rate  equa¬ 
tion  and  stores  the  data  for  output. 

FMDIN 

This  program  sets  the  values  of  the  fracture  mechanics  parameters  AC,  AM,  AP, 
and  AKC  either  from  input  or  from  storage.  The  stored  data  is  set  by  MATTAB 
and  is  contained  in  FMDAT  of  common  block  MATT1. 

FMOUT 

This  subroutine  outputs  the  results  of  the  fatigue,  flaw  growth,  and  residual 
strength  analyses. 

FMTAB 

This  BLOCK  DATA  subroutine  stores  flaw  growth  material  properties  used  in 
the  Erdogan  flaw  growth  equation. 

FRAME 

This  program  sizes  a  frame  based  on  Shanley  criteria  and  minimum  gage  constraints. 
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FUN 


This  subroutine  produces  the  gradient  of  the  objective  function  being  minimized 
at  a  given  design  point  for  use  by  the  nonlinear  math  programming  subroutine 
MINI. 

GEOSTA 

This  program  computes  specific  station  geometry. 

GETMS 

This  subroutine  computes  the  margins  of  safety  for  static  loads  for  stiffened 
panels. 

GETN 

This  subroutine  uses  a  linear- biquadratic  interpolation  scheme  to  find  cycles  to 
failure  given  applied  stress  level  and  stress  ratio  -  SIGMIN/SIGMAX  from  a  set 
of  S-N  data  in  the  form  of  constant  life  diagrams. 

GETS 

This  subroutine  uses  a  linear- biquadratic  interpolation  scheme  to  find  applied 
stress  level  given  stress  ratio  -  SIGMIN/SIGMAX  and  cycles  to  failure  from  a 
set  of  S-N  data  in  the  form  of  constant  life  diagrams. 

GETSTA 

This  subroutine  finds  the  next  station  to  be  optimized  based  on  input  information. 
GINPT1 

This  program  is  called  by  the  input  control  subroutine  INCON.  It  is  used  to  read 
in  the  fuselage  or  aerodynamic  surface  basic  geometry. 

GINPT2 

GINPT3 

GINPT4 

These  are  dummy  programs  which  are  reserved  for  future  alternate  geometry 
input  schemes. 

GROCON 

This  subroutine  determines  the  constant  and  1  G  alternating  stresses  for  each 
segment  of  the  spectrum.  It  also  calls  the  equivalent  spectrum  stress  routine 
and  the  flaw  growth  integration  routine. 
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GSIDE 


This  subroutine  determines  the  constraint  function  for  manufacturing  constraints 
such  as  minimum  gage  and  maximum  stiffener  height,  etc. 

HEADER 

This  program  prints  the  A  PAS- IV  fanfare. 

HSIG 

This  subroutine  computes  the  hoop  stress  distribution  due  to  internal  pressure 
between  adjacent  frames. 

IFRM 

This  subroutine  computes  the  section  properties  of  a  typical  transport  channel  or 
I-section  frame.  A  tear-stopper  inside  die  skin  is  included  with  the  frame  for 
calculating  section  properties. 

INCON 

This  is  an  input  control  program.  This  subroutine  calls  routines  as  required  to 
read  input  data. 

INDEX 

This  subroutine  determines  the  indices  of  the  active  design  variables  for  a 
structural  element. 

INSTIP 

This  subroutine  converts  a  set  of  design  variables  into  a  set  of  detail  geometry 
dimensions.  This  routine  acts  as  an  interpreter  between  the  math  programming 
routine  and  the  structural  analysis  routine. 

INTERP 

This  subroutine  performs  linear  interpolation. 

LDLN1 


This  program  is  called  by  the  input  control  subroutine  INCON.  It  is  used  to 
read  the  fuselage  or  aerodynamic  surface  external  loads. 


LDLN2 

LDLN3 

LDLN4 


These  are  dummy  programs  which  are  reserved  for  future  alternate  external 
load  input  schemes. 

LIFE 

This  program  reads  the  design  criteria  for  fatigue,  flaw  growth,  and  damage 
tolerance  requirements. 


LINK 

This  subroutine  calls  an  input  routine  which  reads  the  information  needed  to  set 
up  the  element  symmetry  groups.  It  then  creates  arrays  of  symmetry  indices. 

LINKED 

This  subroutine  finds  all  of  the  members  of  a  given  symmetry  group. 

LINKIN 

This  subroutine  reads  symmetry  group  information  used  by  LINK. 

LOADS 

This  program  determines  the  externally  applied  loads  at  a  given  station  by  linear 
interpolation  of  input  loads. 

LOADZ 

This  subroutine  determines  the  externally  applied  loads  at  a  given  station  by  linear 
interpolation  of  input  loads. 

LOCALD 

This  subroutine  calculates  the  load  intensities  and  shear  flows  applied  to  the 
structural  elements.  These  applied  internal  loads  are  based  on  the  results  of 
subroutine  BOXLDS  and  the  applied  external  loads  from  subroutine  LOADS. 

LOCATE 

This  subroutine  locates  the  position  of  a  given  station  within  the  stored  geometry 
data  array.  It  then  determines  the  required  interpolation  parameters  needed  to 
extract  the  station  geometry  information. 


LOCBUK 


This  subroutine  evaluates  the  critical  buckling  strain  for  the  stiffened  panel 
s  ub-  elements . 

LOCOPT 

This  subroutine  controls  the  local  design  of  each  symmetry  group.  It  sets  up  the 
input  required  by  the  math  programming  subroutines  and  interprets  the  results. 

MARGIN 

This  subroutine  calculates  the  margins  of  safety  of  composite  panels  using  an 
ultimate  fiber  strain  criteria. 

MATIN 

This  program  reads  the  material  property  input  data. 

MATTAB 

This  BLOCK  DATA  subroutine  stores  the  library  of  material  properties. 

MINI 

This  subroutine  modifies  the  structural  design  of  an  element  to  maximize  the 
margin  of  safety.  The  method  of  Davidon- Fletcher- Powell  is  used. 

MQDL1 

This  subroutine  performs  the  stress  analysis  of  spar-caps  and  longerons. 

NEWTE 

This  subroutine  predicts  the  cross-sectional  area  of  a  structural  element  necessary 
to  produce  a  zero  margin  of  safety. 

NISHEL 

This  program  initializes  the  variables  used  by  the  internal  loads  analysis  routines. 
OFUN 

This  subroutine  is  called  from  ONED  during  the  one  dimensional  search.  It  calls 
FUN  to  evaluate  various  designs. 
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ONED 


This  subroutine  is  called  from  MINI.  It  obtains  the  interval  in  which  a  minimum 
lies  and  performs  a  one-dimensional  minimization. 


OPTCON 

This  program  is  the  optimization  control  routine.  It  does  an  analysis  of  one  struc¬ 
tural  element  or  symmetry  group  at  a  time  until  all  elements  at  a  given  cross- 
section  are  optimized. 

PANPRP 

This  subroutine  computes  the  material  properties  of  a  layered  composite  laminate. 
POEDAT 

This  subroutine  interpolates  from  the  data  stored  in  POET1  to  produce  the  arrays 
needed  by  CRKDAT. 

POETB1 

This  BLOCK  DATA  subroutine  stores  stress  intensity  factor  correction  factors  and 
stiffener  load  concentration  factors  for  cracked  panels  with  riveted  stiffeners. 

PREGRO 

This  subroutine  performs  the  cycle-by-cycle  crack  growth  analysis  and  controls  the 
crack  life  analysis  process. 

PRODAM 

This  subroutine  does  a  Miner's  rule  cumulative  damage  analysis  and  returns  the 
total  damage  and  a  margin  of  safety  in  fatigue. 

PROFLT 

This  subroutine,  based  on  user  directives,  either  extracts  a  flight  spectrum  from 
PROTAB  or  PROTB2,  reads  modification  data  to  these  spectra,  or  reads  a  com¬ 
plete  new  spectrum  from  the  input  data  file. 

PROTAB 

This  subroutine  stores  a  typical  fatigue  spectrum  for  transport  aircraft. 

PROTB2 

This  subroutine  stores  a  typical  fatigue  spectrum  for  a  light-weight  air-to-air 
fighter. 
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PRPMAT 


This  program  defines  the  material  properties  used  during  the  analysis  for  each 
of  the  applied  loading  conditions  at  the  temperature  indicated  for  the  loading 
condition. 


QUAFIT 

This  subroutine  takes  three  values  from  vector  of  X  and  vector  of  Y  starting  at 
first  coordination  point  and  fits  a  quadratic  equation  through  the  three  points  and 
returns  the  value  Y  evaluated  from  the  equation  X. 


RANFIL 

This  subroutine  manipulates  a  random  access  disk  file.  It  may  be  used  to  open, 
read,  or  write  on  the  file. 


RDCMR 

This  is  a  dummy  subroutine  which  may  be  replaced  or  removed  from  the  program 
if  it  is  available  at  the  user's  installation. 

If  the  program  is  available  at  the  user's  installation,  the  program  will  read  the 
control  point  area  for  the  job  and  enable  the  subroutine  TIMIT  to  keep  track  of 
peripheral  processor  time  and  monitor  requests. 

RECFAT 

This  program  controls  the  station  cross-section  redesign  process.  This  process 
attempts  to  minimize  the  station  cross-sectional  area  by  adjusting  the  thickness 
variables  so  as  to  meet  the  constraints. 

RECGRO 

This  program  controls  the  flaw  growth  analysis  and  resizing  procedure  when  flaw 
growth  is  critical.  It  calls  REDSON  to  do  the  analysis  and  NEWTE  to  do  the 
resizing. 

RECOVR 

This  is  a  dummy  subroutine  which  may  be  replaced  or  removed  from  the  program 
if  it  is  available  at  the  user's  installation. 

If  the  program  is  available  at  the  user's  installation,  the  program  will  allow 
recovery  in  the  event  of  a  CP  time  limit,  and  subroutine  CPLIM  will  be  executed. 


RECRES 

This  program  controls  the  residual  strength  analysis  and  resizing  procedure  when 
residual  strength  is  critical.  It  calls  REDSON  to  do  the  analysis  and  NEWTE  to  do 
the  resizing. 
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RE  DC  ON 


This  program  controls  the  resizing  process.  It  calls  the  subroutines  which  add  or 
subtract  material  from  the  structural  elements  in  an  attempt  to  produce  a  minimum 
weight  structure. 

REDOPT 

This  program  produces  an  optimized  cross-section  at  a  station. 


REDSON 

This  subroutine  finds  the  critical  margin  of  safety  for  each  symmetry  group  and 
stores  them. 

RESID 

This  subroutine  determines  the  residual  strength  of  a  panel  for  a  given  crack  size 
and  number  of  broken  stiffeners. 

RIB 

This  program  synthesizes  an  aerodynamic  surface  rib. 

RLIN 

This  is  a  multiple  purpose  linear  interpolation  routine. 

RUNWT 

This  subroutine  calculates  the  running  weights  for  panels,  webs,  and  sparcaps  or 
longerons  at  a  station. 

SAVBT 

This  subroutine  saves  the  B  variables  and  T  variables  to  be  used  as  end  points 
in  the  summary  output  interpolation. 

SEC  MOD 

This  subroutine  calculates  the  secant  modulus  at  a  point  on  the  stress-strain  curve, 
given  the  strain  and  the  three  Ramberg-Osgood  parameters. 

SECPRP 

This  subroutine  computes  the  section  properties  for  a  one,  two  or  three  cell  box 
beam. 
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SETPRO 


This  subroutine  transfers  the  material  properties  of  advanced  composite  materials 
into  the  local  analysis  variables. 

SIGBAR 

This  subroutine  determines  the  equivalent  stress  "SBAR"  and  the  maximum  stress 
"SMAX"  for  a  given  stress  spectrum. 

SKNSTF 

This  subroutine  computes  the  compression  strength  of  dual  material  stiffened 
panels. 

SMINTP 

This  subroutine  interpolates  for  stations  which  were  not  synthesized.  Results  for 
the  interpolated  stations  are  printed  by  subroutine  SUMOUT 

SORT 

This  routine  sorts  the  elements  in  a  given  vector  in  ascending  order.  This  is  a 
COMPASS  routine. 

SPECLD 

This  subroutine  calculates  the  load  intensities  and  shear  flows  applied  to  the 
structural  elements.  These  applied  internal  loads  are  based  on  the  results  of 
subroutine  BOXLDS  and  the  applied  external  fatigue  spectrum. 

STAGE 

This  subroutine  sets  up  station  geometry  when  a  new  station  is  to  be  sized. 
STAOUT 

This  program  prints  results  at  the  end  of  each  station  optimization. 

STATION 

This  program  retrieves  data  from  the  data  bank  for  the  desired  station. 

STORE 

This  subroutine  stores  station  output  information  and  interpolates  for  stations  not 
sized.  It  also  calculates  weights  for  panels,  interior  webs  and  spar-caps  or 
longerons. 
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SUBIN1 


This  subroutine  is  called  by  the  geometry  input  program  GINPT1  to  read  symmetry 
group  input  data. 

SUBIN2 

This  subroutine  is  called  by  the  geometry  input  program  GINPT1  to  read  rib/frame 
input  data. 

SUMOUT 

This  program  is  not  currently  in  use.  It  is  to  be  used  to  output  the  summary 
information  collected  by  STORE. 

TCON 

This  subroutine  initializes  the  geometry  variables  used  by  the  analysis  routines. 
TIMIT 

This  subroutine  logs  the  central  processor  time,  peripheral  processor  time, 
monitor  requests,  and  subroutine  calls  for  multiple  calls  to  multiple  subroutines. 

If  subroutine  RDCMR  is  not  available,  TIMIT  will  not  keep  track  of  peripheral 
processor  time  or  monitor  requests. 

TRP2 

This  subroutine  interpolates  from  the  data  stored  in  BF  to  produce  backface 
correction  factors  for  part  through  cracks  {currently  inoperative). 

WEB1 

This  subroutine  translates  the  optimization  variables  into  detail  geometry  variables 
and  vice-versa.  It  is  used  for  internal  web  elements. 

WEB2 

This  subroutine  performs  stress  analysis  on  internal  shear  web  elements. 

XTOV 

This  subroutine  is  called  by  OPTCON  (entry  VTOX)  to  convert  design  variables  V 
into  optimization  variables  X.  It  is  also  called  by  EVA  to  convert  optimization 
variables  supplied  by  the  math  programming  routine  into  design  variables. 
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SECTION  rv 


USAGE  INSTRUCTIONS 


4. 1  PREPARATION  OF  INPUTS 

4. 1. 1  FORMAT  AND  CONTENT.  The  input  variables  to  the  APAS  program  are 
entered  using  the  NAMELIST  utility.  A  complete  list  of  the  required  inputs  is  pre¬ 
sented  in  the  following  section. 

4. 1. 2  DESCRIPTION  OF  INPUT  VARIABLES.  The  input  data  for  the  APAS  routine 
is  defined  in  Table  4-1.  The  charts  identify  the  variable  names,  descriptions,  and 
units.  The  variables  are  grouped  under  their  appropriate  NAMELIST  input  headings. 

TITLE  (2  LINES) 

$INCN 

$MATIN  (1  SET  FOR  EACH  USER- DEFINED 

MATERIAL) 

$LIF 

$SPEC  (OPTIONAL,  INPUT  ONLY  IF  IN  $UF, 

ID P ROC  -  1  OR  2) 

$FMDM 

$GINPT 

$PARTS 

$LNK 

$LINKN  (1  SET  FOR  EACH  CATEGORY;  I.  E. , 

PANELS,  INTERIOR  WEBS,  STIFFENERS, 
OR  LONGERONS) 

$SUBIN1  (1  SET  FOR  EACH  CATEGORY;  I.  E. , 

PANELS  INTERIOR  WEBS,  STIFFENERS, 
OR  LONGERONS) 

$SUBN2 

ALPHANUMERIC  (1  FOR  EACH  CONDITION.  MUST  HAVE 

LOAD  TITLES  12  LINES) 


$LDLN 


INPUT  DATA  CARDS 
Description:  Title  Cards  (2  cards) 

Format:  8A10 

Column:  1  80 

TITLE  (1)  through  TITLE  (8) 


1 _ 80 

TITLE  (9)  through  TITLE  (16) 


Field  Contents 

TTTLE(I),I»1,  16  Any  alphanumeric  Information  which  the  user  desires  to  input 
for  problem  identification 
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input  data  card 

Description:  User  Supplied  Metallic  Material  Title  Card. 

Format: 

Column: 

Field  Contents 

TYPMAT(I),  Alphanumeric  information  describing  a  metallic  material  to  be 
lxl>g  input  by  the  user. 


8A10 


Wamnrlta; 

1.  These  cards  will  be  repeated  once  for  each  user  defined  metallic 
material. 


INPUT  DATA  CARD 


Description: 

Format: 

Column 

Field 

TTTCON 

NSTAY 

PRESS 


Load  Condition  Header  Card 
4A10,  IS,  5X,  F10. 0 


1  40  41  45  51  60 


TITCON(l)  .  .  .  TTTCON(4) 

NSTAY 

PRESS 

Contents 


Condition  title,  up  to  40  characters 

Number  of  stations  at  which  the  load  case  is  defined  (20  max) 
Internal  pressure  used  for  fuselage  structures,  a  psi 
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TABLE  4-1.  INPUT  DATA  FOR 
PROGRAM:  APAS  IV 


NAMELIST  NAME: 


FILE  NAME: 


_ (Titles  and  Alphanumeric  Information)  _ _ 

SYMBOL  1  DESCRIPTION  SOURCE 

TITLE  (16)  Alphanumeric  Information  for  Problem  Identification*  UNIQ 

TYPMAT 

(10,20)  Alphanumeric  Information  Describing  a  Material  to  be  Input 
3y  User** 

TITCON 

(4,12)  Load  Condition  Title**  UNIQ 


'Located  at  beginning  of  input  deck. 

'*  Located  before  appropriate  input  NAMELISTS 
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TABLE  4-2.  INPUT  DATA  FOR 


PROGRAM:  A  PAS  IV 


NAME  LIST  NAME:  $INCN _  FILE  NAME:  APAS  1 

(Iteration,  Tolerance,  and  Case  Control) _ 


SYMBOL 

DESCRIPTION 

SOURCE 

rri 

Iteration  Count  Limit  on  overall  Redesign/Optimization 

UNIQ 

Procedure  (Default  Value  is  5) 

IT2 

Iteration  Count  Limit  on  fully  Stressed  Redesign  Process 
(Default  is  20) 

IT3 

Iteration  Count  Limit  on  Fletcher  Powell  Optimization 
Procedure  (Default  is  5) 

IT4 

Iteration  Count  Limit  for  Fatigue  Flaw  Growth,  and  Residual 
Strength  Redesign  Process  (Default  is  5) 

EPS1 

Tolerance  on  Redesign  Margins  of  Safety  for  each  Redesign 
Cycle  (Default  is  0.001) 

EPS2 

Tolerance  on  Final  Margins  of  Safety  at  least  one  Non- 
minimum  Gage  Element  of  each  Symmetry  group  for  at  least 
one  load  condition  will  have  a  margin  of  safety  which  satisfies 
MS  EPS2  (Default  is  0.01) 

EPS3 

Tolerance  on  Optimization  Function  Decrease  in  Fletcher 
Powell  Minimization  Technique  (Default  is  .  001) 

EPS4 

Tolerance  on  Design  Variable  Variation  in  Fletcher  Powell 
Minimization  Technique  (Default  is  0.001) 

KEY1 

Specifies  one  of  Four  Geometry  input  subroutines 
=  1,  specifies  general  input  subroutine  "GINPT1" 

=2,3,4  -  not  currently  available 

KEY2 

Specifies  one  of  four  external  loads  definition  subroutines 
=  1  specifies  subroutine  "LDLN1" 

=2,3,4-  not  currently  available 

KEY3 

Specifies  synthesis  of  wing-like  or  fuselage  type  structures 
=  0,  for  fuselage  type  structures 
=>  1,  for  wing-like  structures 

KEY4 

Specifies  number  of  locations  along  structure  to  be  synthe¬ 
sized 

*  1,  synthesize  at  every  rib/frame 

=  2,  synthesize  at  every  second  rib/frame,  etc. ,( maximum 
value  is  50) 

KEY5 

Analysis  /Optimization  flag 
■  0,  perform  optimization 

*  1,  perform  analysis  of  structure  as  input,  do  not  perform 
optimization 

KEY6 

Input  check  flag 

■  0,  normal  mode 

■  1,  check  input  and  quit 

U 

NIQ 

TABLE  4-2.  INPUT  DATA  FOR 
PROGRAM:  APAS  IV 


NAMELIST  NAME:  $INCN  (Cont.) 


FILE  NAME:  APAS  1 


SYMBOL  _ DESCRIPTION _ SOURCE  I 

KEY7  Fatigue  analysis  flag  UNIQ 

=  0,  no  fatigue  analysis 

*»  1,  include  redesign  to  respect  fatigue  criteria 
KEY8  Flaw  growth  flag 

■  0,  no  flaw  growth  analysis 

=  1,  Include  redesign  to  respect  flaw  growth  criteria 
KEY9  Residual  strength  flag 

3  o,  no  residual  strength  analysis 

-  1,  Include  redesign  to  respect  residual  strength  criteria 
KEY10  CP  time  limit  recovery  procedure  flag 

=  0,  use  normal  exist  procedure 
1,  print  output  values  at  time  of  time  limit 
KEY11  Objective  function  derivative  specification  flag  UNIQ 

=  0,  central  difference  derivatives 
»  1,  one  sided  derivatives  (recommended  value) 


TABLE  4-2.  INPUT  DATA  FOR 

PROGRAM:  APAS  IV _ 

NAMELIST  NAME:  S MATIN _ FILE  NAME:  APAS  I 


SYMBOL 

DESCRIPTION 

SOURCE 

NMAT 

MATID(6) 

(Material  specification.  Allows  the  User  to  select  one  of 
the  materials  built  into  the  program  or  to  specify  a  material 
of  his  own. ) 

Number  of  materials  to  be  used 

Identification  number  for  materials  to  be  used: 

ID  Material  Type 

1-3  User  defined  metallic  material 

4  AL-2024-T62 

5  AL-2024-T851 

6  AL-7075-T6 

7  AL-2219-T87 

8  TI-SAL-4V 

9  TI-8AL-1MO-1V 

10  Inconel  718 

11  Inconel  625 

12  RENE'  41 

13  -  IS  User  defined  composite 

16  NARMCO  5505 

17  NARMCO  5206 

\ 

UNIQ 

UNIQ 
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TABLE  4-3.  INPUT  DATA  FOR 

PROGRAM:  APAS  IV _ 

NAME  LIST  NAME:  3  MATIN' _  PILE  NAME:  APAST 


i  VARIABLE 

DESCRIPTION' 

Units 

Data 

FTEN 

Knockdown  factor  applied  to  metallic  material  FTU  to  obtain 

dimen- 

the  allowable  tensile  stress  for  limit  load  conditions 

sionless 

0.66 

NTEMP 

Number  of  temperatures  for  which  a  user  defined  metallic 
material  will  have  temperature  factors  input 

integer 

10 

FTU 

Ultimate  tensile  strength 

psi 

60000.0 

EC 

Modulus  of  elasticity  in  compression 

psi 

10.  5E6 

FCY 

Compressive  yield  strength 

psi 

40000. 0 

FSU 

Ultimate  shear  strength 

psi 

30000. 0 

E 

Modulus  of  elasticity  in  tension 

psi 

10.  3E5 

G 

Shear  modulue 

psi 

3.  3E6 

RHO 

Density 

lb/inJ 

0.  10 

F07 

Stress  from  intersection  of  stress-strain  curve  with  a  secant 
of  slope  0 . 7E  a  Ramberg-Osgood  parameter 

psi 

53000. 0 

EN 

Ramberg-Osgood  shape  parameter  for  the  stress-strain 
curve  -  dimensionless 

dimen¬ 

sionless 

18.  5 

TEMPM(9) 

Temperature  at  which  material  properties  are  being  defined 

*F 

200.0 

FFTU(9) 

Factor  applied  to  the  FTU  at  room  temperature  to  obtain  the 
FTU  at  the  corresponding  temperature  TEMPM(I) 

dimen  - 
sionless 

0.75 

FFCY(9) 

Factor  applied  to  FCY  at  room  temperature  to  obtain  FCY 
at  corresponding  temperature  TEMPM(I) 

dimen¬ 

sionless 

0.75 

FFSU  (9) 

Factor  applied  to  FSU  at.  room  temperature  to  obtain  FSU  at 

dimen- 

0.75 

corresponding  temperature  TEMPM(I) 

sionless 

FEC(9) 

Factor  applied  to  EC  at  room  temperature  to  obtain  EC  at 
corresponding  temperature  TEMPM(I) 

dimen  - 

sionless 

0  -  *  5 

FE(9) 

Factor  applied  to  E  at  room  temperature  to  obtain  E  at 

dimen  - 

0,75 

corresponding  temperature  TEMPM  (I) 

sionless 

FG(9) 

Factor  applied  to  G  at  room  temperature  to  obtain  G  at 
corresponding  temperature  TEMPM(I) 

dimen¬ 

sionless 

0.75 

FRHO(9) 

Factor  applied  to  RHO  at  room  temperature  to  obtain  RHO 
at  corresponding  temperature  TEMPM(I) 

dimen¬ 

sionless 

0.95 

A-I2 


TABLE  4-3.  INPUT  DATA  FOn 

PROGRAM:  A  PAS  IV _ 

NAMELIST  NAME:  SMATIN _  FILE  NAME:  APASI 


VARIABLE 

Ell 

E22 

G12 

U12 

DEN 

EPSAL1 

EPSAL2 

EPSAL3 

EPSAL4 

EPSAL5 


_ DESCRIPTION _ 

Lamina  modulus  of  elasticity  in  the  fiber  direction 
Lamina  transverse  modulus  of  elasticity 
Lamina  in-plane  shear  modulus 

Lamina  Poisson's  ratio  for  loading  in  the  fiber  direction 
Density  of  composite  material 

I  Lamina  ultimate  allowable  tensile  strain  in  the  fiber 
direction 

Lamina  ultimate  allowable  tensile  strain  in  the  transverse 
direction 

Lamina  ultimate  allowable  shear  strain 

Lamina  ultimate  allowable  compressive  strain  in  the  fiber 

direction 

Lamina  ultimate  allowable  compressive  strain  in  the 
transverse  direction 


Units 


Data 


psi 

psi 

psi 

dimen¬ 

sionless 


40.0E6 

1.5E6 

4. 0E6 

0.20 


lb/in3 


0.  55 


in. /in.  7.0E-3 

in. /in.  7.0E-3 


rad.  r 

27. 0E  -3 

in. /in. 

11. 0E  -3 

in. /in. 

28.0E-3 
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TABLE  4-4.  INPUT  DATA  FOR 


PROGRAM: 


APAS  IV 


NAME  LIST  NAME: 


FILE  NAME: 


APASI 


(Spectrum  Input  Data) 


VARIABLE 

DESCRIPTION 

UNITS 

(Spectrum  modification  is  required  if 

ID  PROC  =  1  or  2  in  NAMELIST  $LIF.  If 

ID  PROC  =  2,  the  following  variables  would 
describe  the  fatigue  spectrum.  One  load 
condition  in  NAMELIST  $LDLN  would  define 
the  reference  fatigue  leading  condition. ) 

NOCYC 

Number  of  load  steps  in  the  input  stress 
spectrum 

Integer 

EM 

Number  of  flights  represented  in  the  input 
stress  spectrum 

Real 

FSIG 

Factor  applied  to  FKMAX  and  FKMIN 

Real 

FKMAX(IOOO) 

Factor  applied  to  reference  fatigue  condition 
stress  to  obtain  maximum  spectrum  stress 
for  each  of  NOCYC  load  steps.  Maximum 
spectrum  stress  FMAX,  based  on  reference 
stress  SIG,  is  obtained  as  follows: 

Real 

FMAX  =  FKMAX  *  FSIG  *  SIG 

FKMIN(IOOO) 

Factor  applied  to  reference  fatigue  condition 
stress  to  obtain  minimum  spectrum  3tress 
for  each  NOCYC  load  steps 

Real 

CYC  (1000) 

Number  of  load  cycles  for  each  of  NOCYC  load 
steps 

Real 

L5 


TABLE  4-4  INPUT  DATA  FOR 


PROGRAM:  APAS  IV 


NAME  LIST  NAME:  SSPEC _  FILE  NAME:  APASI 


(Spectrum  Input  Data) 


VARIABLE 

DESCRIPTION 

UNITS 

DATA 

CY(20,20) 

Cycles  CY(J,  I)  in  the  profile  for  segment  I, 
subsegment  I  based  on  CYCBS  Bights 

Real 

FCON(6,20) 

Constant  stress  composition  table,  FCON(L,  I), 
where  the  constant  stresses  for  segment  I  are 
based  on  linear  combinations  of  the  stresses 
due  to  L-spectrum  loading  conditions.  Spec¬ 
trum  loading  conditions  are: 

Real 

Cond  No.  Description 

1 

1  1  G  taxi 

2  1  G  Bight 

3  1  G  flight  + 1  G  vertical  gust 

4  1  G  Bight  -  1  G  maneuver 

5  1  G  landing  impact 

6  Maximum  internal  pressure 


FALT(6,20) 


I 


Alternating  stress  composition  table, 

FALT(L,  I),  where  stress  excursions  for  seg¬ 
ment  I  are  based  on  the  linear  combination  of 
the  stresses  due  to  L-spectrum  loading  condi¬ 
tions-  Magnitude  of  the  stress  excursions 
is  obtained  from  linear  scaling  by  the  incre¬ 
mental  load  factor  table  DG,  based  on  the 
load-type  specification  LT. 


Real 


j 

i 


r 


TABLE  4-4  INPUT  DATA  FOR 

PROGRAM:  A  PAS  TV _ 

NAME  LIST  NAME:  SSPEC _  FILE  NAME:  APASI 


(Spectrum  Input  Data) 


VARIABLE 

DESCRIPTION 

UNITS 

DATA 

(Spectrum  modification  is  required  if 

ID  PROC  =  1  or  2  in  NAMELIST  $LIF.  If 

ID  PROC  =  1  the  following  variables  may  be 
used  to  modify  stored  data. ) 

CYCBS 

Basis  for  occurrence  data,  based  on  CYCBS 
flights 

Flights 

10, 000 

NOSEG 

Number  of  segments 

Integer 

20 

NOSUB 

Number  of  subsegments 

Integer 

20 

DG(20) 

Incremental  load  factor  for  each  subsegment 

Real 

LT(20) 

Segment  load- type  specification  describing 
nature  of  cycling 

Integer 

1  =  maneuver;  twice  delta  G  excursion 
from  basic  condition 

2  =  taxi/gust;  positive  and  negative 
delta  G  excursion  about  basic 
condition 

! 

i 

! 

3  =  landing  impact;  twice  delta  G 
negative  excursion  from  basic 
condition 

i 

IN (20) 

Segment  environment  specification 

Integer 

i 

1  =  flight  condition 

i 

i 

2  =  ground  condition 

4-17 


TABLE  4-5  INPUT  DATA  FOR 


1 


PROGRAM:  APAS  IV _ 

NAME  UST  NAME:  SFMDN _  FILE  NAME:  _ A  PA  SI 


(Fracture  Mechanics  Material  Data) 


VARIABLE 

DESCRIPTION 

UNITS 

DATA 

IFM(6) 

Flag  to  indicate  whether  to  use  stored  fracture 
mechanics  material  properties  data  or  to 
input  the  data 

Integer 

0  =  stored  data  will  be  used 

1  =  user  inputs  data 

AC  (6) 

Growth- rate-equation  coefficient  (c) 

Real 

1.0E-20 

AM(6) 

Growth-rate-equation  (1-R)  exponent  (m) 

Real 

0.6 

AN(6) 

Growth-rate-equation  exponent  (n) 

Real 

3.64 

AQ<6) 

Acceleration  model  exponent  (q) 

Real 

0.3 

AKC(6) 

Plane  stress  facture  toughness 

psi  s[ In. 

9.2E4  1 

AKIC(6) 

Plane  strain  fracture  toughness 

psi  sT in. 

4.5E4  | 

1 

THRESH(6) 

Threshold  value  of  Ak  at  R  =  0 

psi  \T in. 

2.5E3  1 

CUT(6) 

Shutoff  ratio 

Real 

1 

3.0  i 

| 

RCUT(6) 

Positive  stress  ratio  cut  off  value 

Real 

0.75  | 

i 

RCUTN(6) 

Negative  stress-ratio  cutoff 

Real 

i 

-0.99  j 

i 

I 

j 

I 

1 

1 

! 

i 

4-13 


TABLE  4-6.  INPUT  DATA  FOR 


NAMELIST  NAME: 


PROGRAM: 

SGINPT 


APAS  IV 


FILE  NAME:  APASI 


(Box  Beam  Geometry  Input) 
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VARIABLE 

DESCRIPTION 

Units 

Data  1 

NODES 

Number  of  node  points  in  a  cross  section. 

3  *  nodes  *  20 

integer 

is  : 

I 

NWEB 

Number  of  interior  webs  maximum  of  3 

integer 

2  | 

NLONG 

Number  of  spar  caps /longerons  max  of  10 

integer 

5  1 

O  1 

NSTAG 

Number  of  geometry  control  stations  max  of  20 

integer 

15  I 

STAG{20) 

Station  number  of  control  station 

in. 

FRSP(20) 

Rib/Frame  spacing 

in. 

XLDRF(20) 

X-Coordinate  of  input  loads  reference  axis 

in. 

200.0 

ZLDRF(20) 

Z-Coordinate  of  input  loads  reference  axis 

in. 

GTRX(20) 

Taper  distance  of  structure  in  the  X- direction 

in. 

100.0 

GTRZ(20) 

Taper  distance  of  structure  in  the  Z- direction 

in. 

ITEM  (20) 

Node  number.  Begins  with  1,  numbered  clockwise 
around  structure  cross-section 

integer 

15 

GX(20,  20) 

GX(I,  J)  =  X  Coordinate  of  Node  J  at  Station  I 

in. 

100.0 

GZ(20,  20) 

GZ(I,  J)  =  Z  Coordinate  of  Node  J  at  Station  I 

in. 

rw  (6 ) 

IW(I)  =  Number  of  the  first  node  to  which  the  first 
interior  web  is  attached.  IW{6)  =  number  of  the  node 
to  which  the  other  end  of  the  first  interior  web  is 
attached.  IW{2)  and  IW(5)  describe  the  next  interior 
web,  IW(3)  and  IW{4)  describe  the  last  interior  web. 

integer 

2  | 

1 

IL(10) 

Node  number  of  spar-cap/longeron  I 

integer 

5  ; 

DL(10) 

Orientation  angle  of  spar-cap/longeron  I 

degrees 

20.  0 

1 

i 

± — A- 


J 


THEORETICAL 

intersection  geometry 


VIEW  LOOKING  AFT 


outboard 


geometry 


Figure  4-3.  Typical  Wing  Taper  Ratio 
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TABLE  4-7.  INPUT  DATA  FOR 


PROGRAM:  APAS  IV 


NAMELIST  NAME :  SLNK _  FILE  NAME:  APASI 

(Symmetry  Group  Control) 


VARIABLE 

DESCRIPTION 

Units 

Data 

NSGP 

Number  of  panel  symmetry  groups 

integer 

5 

NSGW 

Number  of  interior  web  symmetry  groups 

integer 

5 

NSGL 

Number  of  sparcap/ longeron  symmetry  groups 

integer 

5 
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NSGP 


Number  of  panel  symmetry  groups 


NSGW 

NSGL 

Remarks: 


Number  of  interior  web  symmetry  groups 

Number  of  sparcap/longeron  symmetry  groups 

1.  The  use  of  symmetry  groups  is  illustrated  in 
Figure  4-4. 

2.  The  following  examples  illustrate  the  use  of  symmetry  groups. 

The  first  example  represents  a  section  cut  through  the  fuselage 
of  a  typical  transport  fuselage.  The  numbers  indicate  individual 
panel  elements  between  adjacent  node  points.  A  typical  set  of 
symmetry  groups  is  listed  below  the  fuselage  section,  indicating 
that  four  separate  designs  are  desired,  one  for  each  group,  i.e. , 
panels  1,  2,  3,  16,  17,  and  18  will  all  have  the  same  detailed 
design  dimensions.  This  design  will  be  dictated  by  the  most 
critical  panel  or  panels  in  the  group.  Note  that  corresponding 
panels  on  opposite  sides  of  the  centerline  are  members  of  the 
same  symmetry  group,  (e.  g. ,  panels  1  and  18  are  members  of 
Group  2).  This  type  of  grouping  provides  for  centering  symmetry. 

The  second  example  presents  the  use  of  symmetry  groups  for  a 
wing-type  structure.  As  before,  the  numbers  indicate  panel 
numbers,  with  the  exception  of  Wl,  which  represents  an  interior 
web  member.  The  set  of  symmetry  groups  shown  below  the  wing 
section  indicates  four  symmetry  groups.  Panels  7,  14,  and 
interior  web  Wl  are  not  members  of  any  symmetry  group  for  this 
example  and  hence  will  have  unique  designs.  Therefore,  a  total 
of  seven  separate  designs  will  be  generated  at  each  cross -section. 
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A  new  set  of  designs  for  a  cross-section  is  produced  at  each 
station  where  optimization  is  performed.  The  symmetry 
grouping  is  preserved  for  the  entire  length  of  the  structure. 
Panels  of  a  given  symmetry  group  will  have  the  same  design  at 
any  given  station.  However,  the  design  is  free  to  change  from 
station  to  station. 


wni 

uian.il 

ll.tlU 

titan 


la^iLTnriWiMi 


Figure  4-4.  Structural  Symmetry  Grouping 


TABLE  4-8.  INPUT  DATA  FOR 


PROGRAM:  A  PAS  IV 


NAMELIST  NAME:  $LINKN 


FILE  NAME:  APASI 


VARIABLE 


DESCRIPTION 


KMAX(20)  KMAX(l)  =  Number  of  elements  in  symmetry  Group  I  integer 

NSE(20,  20)  NSE  (K,  I)  a  Element  number  of  the  Kth  element  in  integer 
symmetry  Group  I 

NOTE:  One  L1NKIN  namelist  is  input  for  each 
symmetry  group  type  -  panels,  webs,  and  spar-caps 
/longerons.  Total  of  3. 

Numbering  of  panels  is  clockwise  around  the  cross 
section.  Numbering  of  webs  starts  with  web  closest 
to  front  spar.  Spar-cap/ longeron  numbers  are  equal 
to  their  corresponding  node  numbers. 


& 


TABLE  4-9.  INPUT  DATA  FOR 

PROGRAM:  APAS  IV _ 

NAMELIST  NAME :  SSUBIN1 _  FILE  NAME: _ 4 

(Structural  Element  Configuration  Type  Identification 


VARIABLE  DESCRIPTION 


ITYPE(20)  Structural  configuration  type  number. 

IDSET(20)  Unique  ID  for  the  corresponding  structural  element 
configuration  specification. 

IDMAT1(20  Material  ID  of  element  for  riveted  panel  elements. 
This  is  ID  of  skin  only. 

IDMAT2(20  Material  ID  of  stiffeners  for  riveted  panel  only: 

ITYP(I)  =  4  thru  9. 

FELEM(20  Used  for  panel  and  web  element  types  10,  11,  and  12 
only.  Unsupported  panel  width  in  terms  of  a  factor 
times  the  panel  element  width. 

T(20,  4)  Initial  value  for  various  cross-section  thicknesses 
(T  variables). 

TMN(20,  4)  Minimum  values  for  T  variables. 

B(20,  4)  Initial  values  for  various  cross-section  dimensions 
(B  variables). 

BMN20,  4)  Minimum  values  for  B  variables. 

BMX(20, 4)  Maximum  values  for  B  variables. 


Units 

Data 

integer 

1 

integer 

2 

integer 

2 

integer 

2 

In. 

20.0  ' 

in. 

.04 

in. 

.04 

in. 

10.0 

in. 

5.0 

in. 

100.0 
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TABLE  4-10.  INPUT  DATA  FOR 

PROGRAM:  APAS  IV _ 

NAMELIST  NAME :  SSUBN2 _  FILE  NAME: 

(Rib /Frame  Configuration  Type  Identification) 


SYMBOL  DESCRIPTION 


Rib /Frame  configuration  type  number 
Frame 

=  0,  suppress  frame  analysis 
=  I».  ring  frame  with  zee  cross  section 


A  PA  SI 


SOURCE 


UNIQ 


IDMTFR 

IDMTRS 


3  0,  suppress  rib  analysis 
=  1,  corrugated  web 

-  2,  integral  web 
=  3,  built-up  web 

-  4,  integral  truss 

-  5,  built-up  truss 

Rib/Frame  material  ID 

Rib  stopper  material  ID  used  for  fuselage  type 
structures  only 


UNIQ 

UNIQ 
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TABLE  4-11.  INPUT  DATA  FOR 
PROGRAM:  APAS  IV 
NAMELIST  NAME:  SLDLN  _ 


FILE  NAME: 


A  PAST 


V  A  RUBLE 

DESCRIPTION 

Data 

NCOND 

Number  of  loading  conditions  maximum  of  6 

integer 

2 

FULT 

Ultimate  factor  of  safety 

real 

1. 5 

NSPEC 

Number  of  loading  conditions  in  the  fatigue  and  flaw 
growth  spectrum 

integer 

5 

FBUCK 

Factor  for  initial  buckling  of  skin  between  stiffeners 
for  skin  stiffener  construction  type  panels 

real 

0.  50 

NSTAY{12; 

Number  of  stations  at  which  load  case  is  defined 

integer 

10 

PRES(12) 

Internal  pressure  used  for  fuselage  structures 

psi 

15.0 

ELIN(12 ) 

Factor  applied  to  STA 

real 

0.50 

FLD(12) 

Factor  applied  to  all  input  load  components 
(temperature  excluded) 

real 

1.  5 

FA(12) 

Factor  applied  to  AX 

real 

1.5 

FXS(12) 

Factor  applied  to  XS 

real 

1. 5 

FZS(12) 

Factor  applied  to  ZS 

1.5 

FT0R(12) 

Factor  applied  to  TOR 

1.5 

FXM(12) 

Factor  applied  to  XMOM 

real 

1.  5 

F  ZM(12 ) 

Factor  applied  to  XMOM 

real 

1. 5 

FTEMP(12 

)  Factor  applied  to  TEMP 

real 

1.5 

STA(20,  12 

Fraction  of  fuselage  length  or  wing  semi-span 

real 

0,  50 

AX(20.  12) 

Station  axial  load 

lbs 

10000. 0 

XS(20,  12) 

Station  X  shear  force 

lbs 

10000. 0 

ZS(20,  12) 

Station  Z  shear  force 

lbs 

10000. 0 

TOR(20,  12 

) Torque 

in- lbs 

10000. 0 

XMOM(20. 

12)  Bending  moment  about  X-axis 

in-  lbs 

10000.0 

ZMOM(20, 

12)  Bending  moment  about  Z-axis 

in-  lbs 

10000.0 

TEMP(20, 

•  2)  Structural  temperature 

deg/Fahr 

200.  0 
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Figure  4-5  shows  a  graphical  representation  of  the  flow  of  data  into  and  out  of  the 
A  PAS  Structural  Synthesis  module.  The  NAMELIST  input  data  is  on  a  file  called 
APASI.  This  data  is  retrieved  from  the  data  bank  by  the  ADM  and  preprocessed 
to  the  forma.t  needed  by  the  sldn  panel  module.  The  output  is  stored  on  files  called 
APAS<D1  and  APA502.  This  output  data  file  is  post-processed  by  the  ADM  and  stored 
in  the  data  base.  The  output  is  discussed  in  the  OUTPUT  FORMAT  and  CONTENT 
section  of  this  manual . 


Figure  4-5.  A  PAS  -  Data  Flow 
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4.1.3  LIMITATIONS  AND  RESTRICTIONS 

The  APAS  program  performs  a  multi- station  structural  analysis  for  sizing  of  box 
beam  structure  element  A  detailed  discussion  of  the  analysis  procedure  and 
assumptions  Is  presented  in  the  following  paragraph. 


TECHNICAL  DISCUSSION 


The  technical  approach  used  in  APAS  IV  is  applicable  to  any  closed  section  beam-like 
structure,  and  it  is  typical  of  the  procedure  used  in  the  early  design  phase  of  aircraft 
stiucture.  The  overall  approach  makes  use  of  a  point  design/analysis/redesign 
process  that  is  iterated  until  an  acceptable  design  is  produced.  Figure  2-1  presents 
the  functional  flow  chart  for  the  approach  used  in  the  APAS  IV  computer  program.  This 
flow  chart  outlines  the  major  analysis  and  design  loops  of  the  program. 

This  section  includes  a  discussion  of  each  of  the  following  topics:  Component  Geometry, 
Structural  Elements,  Flight  Profile  ard  Load  Spectrum,  External  Loads,  Structural 
Design  Procedure,  and  Structural  Analysis.  The  structural  analysis  discussion 
includes  the  static  strength,  stability,  fatigue,  fracture,  and  residual  strenght  analysis 
used  in  APAS  IV. 

COMPONENT  GEOMETRY 

The  geometry  of  each  component  (fuselage,  wing,  horizontal  and  vertical  stabilizer)  is 
represented  by  the  coordinates  of  a  set  of  nodes  at  each  of  the  various  stations  along 
the  component.  This  nodal  geometry  describes  the  shape  of  the  component  used  for 
the  computation  of  section  properties  and  internal  loads.  The  program  is  capable  of 
reading  and  storing  nodal  information  at  geometry  control  stations.  Storage  for  20 
control  stations  is  available  for  nodal  information;  however,  fewer  may  be  used.  The 
program  uses  linear  interpolation  between  control  stations  to  determine  required 
nodal  information.  Nodal  information  at  a  control  station  consists  of  X  and  Z 
coordinates  for  each  node.  The  program  provides  for  a  maximum  of  20  nodes  per 
control  station. 

FUSELAGE  NODAL  GEOMETRY 

Nodal  geometry  for  a  typical  transport  fuselage  is  presented  in  Figure  4-6.  Nodes 
are  numbered  starting  at  the  top  centerline  and  proceeding  clockwise  looking  aft. 

AERODYNAMIC  SURFACE  NODAL  GEOMETRY.  The  nodal  geometry  describes 
the  box  structure  for  an  aerodynamic  surface  with  up  to  five  spars.  The  nodes 
are  numbered  beginning  at  the  upper  sparcap  of  the  front  spar  and  proceeding  clock¬ 
wise  to  the  lower  front  sparcap.  A  typical  surface  nodal  geometry  is  presented  in 
Figure  4-7. 


4-30 


Section  A-A  (Typ) 

Figure-  4-6.  Fuselage  Nodal  Geometry 
STRUCTURAL  ELEMENTS,  RIBS,  AND  FRAMES 

STRUCTURAL  ELEMENTS.  Structural  elements  Include  skin  panels,  spar 
webs,  and  spar  caps.  Each  element  is  described  by  a  type  number  and  by  from  one 
to  eight  dimension  variables.  The  dimension  variables  are  of  two  types,  thickness 
variables  and  non-thickness  variables  such  as  stiffener  spacing,  stiffener  height,  and 
corrugation  angle.  All  variables  may  have  manufacturing  constraints  imposed.  In 
general,  non-thickness  variables  (i.  e. ,  B  variables)  may  be  set  to  a  constant  value 
or  may  be  constrained  between  upper  and  lower  limits.  Thickness  variables  (1.  e. , 

T  variables)  may  either  be  set  to  a  constant  value  or  may  have  minimum  gage  con¬ 
straints  imposed.  T  and  B  variables  are  shown  for  the  various  construction  types  in 


Figure  4-7 ;  Aerodynamic  Surface  Nodal  Geometry 


Skin  Panel  Construction  Types.  The  structural  synthesis  program  Includes 
12  types  of  panel  elements  as  presented  in  Figure  4-8  The  stiffeners  on  panel  types 
one  through  nine  are  assumed  to  be  oriented  parallel  to  the  elastic  axis  of  the  struc¬ 
ture.  The  0  degree  ply  of  panel  type  12  is  also  assumed  to  he  parallel  to  the  elastic 
axis. 


"Spar  Web"  Elements.  The  structural  synthesis  program  contains  seven 
types  of  "spar  web"  elements.  Four  of  these  are  truss  type  elements,  two  are  stiffen¬ 
ed  webs,  and  the  remaining  one  is  a  corrugated  web.  These  elements  are  presented  in 
Figure  4-9  .  "Spar  Web"  elements  are  assumed  to  resist  only  shear  and  crushing 
loads,  the  axial  stiffness  of  these  elements  is  assumed  to  be  zero  for  the  purpose  of 
computing  section  properties. 


Figure  4-9.  Spar  Web  Elements 
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Spar  Cap  Elements.  Four  types  of  spar  caps  are  currently  available.  They 
include  integral  tee  and  angle  and  riveted  tee  and  angle  as  shown  in  Figure  4-10. 
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Figure  4-10.  Spar  Cap  Elements 


RIBS.  The  types  of  ribs  available  within  the  program  are  presented  in  Figure  4- 
11.  The  ribs  consist  of  caps  and  webs  or  truss  elements.  Rib  caps  are  sized  to  re¬ 
act  a  moment  at  the  rear  spar  due  to  the  loading  on  the  surface  aft  of  the  rear  spar. 
Rib  webs  are  sized  to  carry  shear  and  to  support  crushing  loads. 


built-up 

WEB 


BUILT-UP 

TRUSS 


CORRUGATED 

WEB 


Figure  4-11 


FRAMES.  A  typical  ring  frame  is  shown  in  Figure  4-12.  The  frames  are 
sized  so  that  the  outer  flange  clears  all  of  the  skin  stiffeners.  The  inner  flange  is 
maintained  at  14  cm  (5. 5  inches)  from  the  outer  skin  contour.  The  frame  is  sized 
using  Shanley's  criteria  to  set  a  minimum  frame  bending  stiffness.  The  frame  is  set 
to  mini  mum  gage  for  non-critical  areas. 


El  « 


Shanley's  criteria  (Reference  4) 


where: 

El  *  frame  bending  stiffness  [in2  -  lb.  ] 

M  ■  maximum  resultant  fuselage  bending  moment,  s/M^2  +  M^2  [in.  -lb.  ] 
D  »  fuselage  diameter  [In.] 

L  **  frame  spacing  [in.  ] 

Cj  *  fit  coefficient  (.00025)  [dimensionless] 


Figure  4-12  Typical  Ring  Frame 


FLIGHT  PROFILE  AND  LOAD  SPECTRUM 


The  fatigue  load  spectrum  defines  the  number  of  times  that  incremental  loads  of  given 
magnitudes  are  encountered  during  the  design  life  of  the  aircraft.  Experimental  data 
is  available  that  defines  the  probable  magnitudes  and  frequency  of  occurrence  of  these 
incremental  loads  as  a  function  of  aircraft  type,  configuration  parameters,  and  flight 
parameters. 

The  configuration  and  flight  parameters  are  defined  using  a  typical  flight  profile,  which 
is  divided  into  segments.  Parameters  are  averaged  for  each  segment,  and  these  aver¬ 
age  values  are  used  in  finding  the  incremental  loads.  See  Figures  4-13  and  4-14. 

Typical  flight  spectra  based  on  medium  range  operation  of  a  contemporary  transport 
aircraft  and  for  a  light  weight  air-to-air  fighter  are  currently  included  in  the  program 
library  for  fatigue  and  flaw  growth  analysis. 

The  parameter  values  for  each  transport  aircraft  flight  segment  are  listed  in  Table 
4-12.  The  segments  are  divided  into  subsegments,  with  each  subsegment  represent¬ 
ing  a  particular  magnitude  of  incremental  load.  Using  the  segment  parameters  and 
the  subsegment  load,  frequency  of  occurrence  of  the  incremental  load  is  found  for  each 
subsegment  using  the  methods  and  information  in  Reference  5. 

For  gust  loads,  curves  showing  gust  velocity  vs  frequency  of  occurrence  are  found  in 
Reference  5,  Figures  C 13-32  through  C 13-37.  From  Reference  5,  Page  C13  24, 


Ag  =  mSVeUdeKg  Po/<2W) 

For  speeds  below  critical  Mach  number: 

K  ,  •88|‘8  „  ,  _2W_ 

g  5.3+Hg  Hg  mgeSp 

where 

Ag  -  incremental  load  factor  [dimensionless] 

m  =  slope  of  lift  curve  [dimensionless] 

2 

S  =  wing  area  [Ft.  ] 

Vg  =  equivalent  airspeed  [Ft.  /Sec.  ] 

Ude  =  derived  gust  velocity  [Ft.  /Sec.  ] 

K  =  gust  alleviation  factor  [dimensionless] 

O  Q 

pQ  =  air  density  at  sea  level  [Ib/ft  ] 

W  =  aircraft  weight  [lb.  ] 

p  =  aircraft  mass  ratio  [dimensionless] 
o  2 

g  =  acceleration  of  gravity  [Ft. /Sec.  ] 

c  =  mean  geometric  wing  chord  [Ft.  ] 

P  =  air  density  [lb.  /ft.  ®] 
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OCCURRENCES  PER  FLIGHT  HOUR 


SEGMENT 

Figure  4-13.  ‘typical  Flight  Profile 


5  10 

INCREMENTAL  LOAD  FACTOR  (A?) 


Figure  4-14.  typical  Segment  Load  Frequency  Curve 
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Solving  for  the  derived  gust  velocity, 


TJ  m 

de  mS  V  K  p  /(2W) 

ego 

Ude  is  then  calculated  for  each  subsegment,  and  the  curves  of  Figure  C13-37  of  Ref¬ 
erence  5  are  used  to  find  the  frequency  of  occurrence. 

For  maneuver  loads.  Figure  C13-41  of  Reference  5  shows  Incremental  load  factor 
versus  frequency  of  occurrence.  For  taxi  loads.  Figure  C13-46  of  Reference  5  shows 
incremental  load  factor  versus  frequency  of  occurrence.  Incremental  load  factor  ver¬ 
sus  frequency  of  occurrence  for  landing  loads  was  averaged  from  data  for  two  com¬ 
mercial  transport  aircraft. 

The  resulting  fatigue  load  spectrum  is  shown  in  Table  4-13.  The  number  of  cycles  is 
based  10, 000  flights.  The  variation  between  cycles  and  flights  is  linear,  so  that 
linear  ratloing  of  cycles  and  design  life  is  valid. 


Table  4-12.  Typical  Transport  Flight  Profile 
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The  flight  profile  used  as  the  basis  for  developing  a  typical  fighter  service  load 
spectrum  is  presented  in  Table  4-14  This  profile  presents  the  segments  of  an  air-to- 
air  combat  mission  for  a  typical  lightweight  fighter.  Maneuver  loads  for  the  flight 
segments  were  obtained  from  Table  I  of  reference  6,  which  presents  the  data  in  the 
form  of  load  factor  versus  cumulative  occurrences.  In  conjunction  with  the  foregoing 
data,  a  representative  supersonic  air-to-air  combat  spectrum  has  been  added  as 
shown  in  Table  4-15.  Taxi  load  factor  versus  cumulative  occurrence  data  was  obtained 
from  Table  Vin  of  reference  6.  Incremental  load  factor  versus  frequency  of  occur¬ 
rences  for  the  landing  segment  was  derived  from  the  sink  speed  versus  landing  sink 
speed  was  converted  to  vehicle  load  factor  by  assuming  a  landing  gear  stroke  of 
12  inches  and  oleo  efficiency  factor  of  0.  8. 


Table  4 -14 

TYPICAL  LIGHTWEIGHT  FIGHTER  AIR-TO-AIR  COMBAT  MISSION 


Segment 

Gross  Weight 
(lb) 

Altitude 
(1,000  ft) 

Mach  No. 

Distance 
(n  mi) 

Time 

(min) 

1. 

Taxi 

12  355 

SL 

- 

- 

- 

2. 

Climb 

12,855  -  12,480 

0-46 

0.  72  -  0.  85 

42.6 

5.364 

3. 

Cruise 

12,480  -  12,175 

46 

0.  85 

112.3 

13. 820 

4. 

Combat 

12,175  -  12,005 

30 

0.90 

- 

0.620 

5. 

Accelerate 

12,005  -  11,740 

30 

0.9  -  1.4 

8.15 

0.707 

6. 

Combat 

11,740  -  10,260 

30 

1.4 

- 

1.125 

7. 

Cruise 

10,260  -  9,920 

50 

0.85 

154.9 

19.290 

8. 

Loiter 

9,920  -  9,625 

10 

0.33 

- 

15.0 

9. 

landing 

9,625 

SL 

- 

- 

- 
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TABLli  4-16  TYPICAL  AIR-TO-AIR  IICUITR  PATIGUI:  SPtiCWLM  -  CYC141S  PIR  1,000  PLUSH'S 


external  loads 

Net  limit  loads  due  to  the  air  loads,  inertia  loads,  and  landing  gear  loads  for  various 
flight  and  ground  conditions  are  input  to  the  program. 

The  loading  conditions  are  separated  into  two  groups.  The  first  group  consists  of 
from  one  to  six  conditions.  These  conditions  are  specified  by  the  user  and  are  used 
to  size  the  structure  so  as  to  preclude  static  strength  failures  and  to  meet  residual 
strength  requirements.  Hie  second  group  consists  of  the  six  conditions  listed  in 
Table  4-17.  These  conditions  are  used  to  define  the  fatigue  stress  spectrum  described 
used  in  the  structural  analysis. 


Table  4-17.  Fatigue  Spectrum  Loading  Conditions. 


Condition 

Number 

i 

Description 

1 

1G  Taxi 

2 

1G  Level  Flight 

3 

2G  Vertical  Gust 

4 

2G  Maneuver 

5 

1G  Landing  Impact 

6 

Maximum  Pressure 

(Fuselage) 

Each  loading  condition  defines  the  six  components  of  load  (AX,  XS,  ZS,  TOR,  XM, 
ZM)  at  up  to  20  stations  along  the  structure.  The  sign  convention  used  is  presented 
in  Figure  4-15.  A  typical  fuselage  loading  condition  is  illustrated  in  Figure  4-16 
Steps  in  the  loading  curves  are  represented  by  repeating  stations  with  the  two  differ¬ 
ent  load  component  values.  The  reference  axis  used  for  input  loads  is  the  centerline 
for  fuselages  and  line  midway  between  the  front  and  rear  spars  for  aerodynamic 


Figure  4-16.  Typical  Fuselage  Load  Condition 


STRUCTURAL  DESIGN  PROCEDURE 

The  structural  design  procedure  starts  with  the  input  design,  then  through  a  series  of 
design  analysis  and  redesign  Iterations  produces  a  final  design, which  satisfies  the  de¬ 
sign  criteria.  The  Iteration  process  continues  until  two  successive  iterations  produce 
a  change  in  weight  that  is  within  a  specified  tolerance. 

The  design  analysis  Involves  the  comparison  of  applied  stresses  and  allowable  stresses. 
The  Internal  load  solution  Is  used  to  calculate  the  applied  stresses.  The  box  beam  lae 
Internal  load  solution  was  selected  Instead  of  a  finite  element  solution  In  order  to 
keep  computer  execution  time  at  an  acceptable  level.  This  selection  restricts  this 
procedure  to  relatively  clean  beam-lOce  structures.  However,  decoupling  of  Internal 
loads  from  one  station  to  the  next  Is  basic  to  the  box  beam  theory.  Hence  the 
overall  design  problem  Is  reduced  to  a  series  of  cross-section  design  problems 
at  discrete  locations  along  the  structure. 
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The  procedure  used  to  design  the  cross-section  is  a  two-part  procedure.  The  first 
part,  the  section  sizing  procedure,  adds  or  subtracts  material  from  the  structural 
elements  of  the  cross-section  in  order  to  produce  a  zero  margin  or  minimum  gage 
design.  The  second  part  employs  a  non-linear  programing  technique  to  maximize  the 
efficiency  of  each  element  while  maintaining  a  constant  weight  design.  This  element 
optimization  procedure  is  then  iterated  with  the  section  sizing  procedure  until  the  de¬ 
sign  converges.  Convergence  occurs  when  two  successive  iterations  produce  a  change 
in  weight  that  is  within  a  specified  tolerance. 

SECTION  SIZING  PROCEDURE.  The  section  sizing  procedure  sizes  the  structure 
based  on  design  criteria  such  as  static  strength,  stability,  service  life,  and  residual 
strength. 

The  procedure  used  to  size  the  structure  is: 

a.  Analyze  the  structure  as  it  is  defined  by  the  input  data. 

b.  Predict  new  skin  thickness  and  stiffener  area  based  on  the  analysis  results. 

c.  Re-analyze  the  structure  as  predicted  is  b. 

Steps  b.  and  c.  are  iterated  until  the  minimum  weight  structure  satisfying  the  design 
criteria  is  found.  During  this  process,  material  is  added  or  removed  from  the  panel 
such  that  the  design  proportions  produced  in  the  optimization  phase  are  maintained. 

The  equivalent  thickness  (t)  of  a  structural  panel  is  computed: 


where:  t^  -  skin  thickness  [In.  ] 

A^  -  stiffener  area  [in.2] 

-  stiffener  spacing  [In.  ] 

-  stiffener  material  density  [lb.  /in. 3] 

Pa^  -  skin  material  density  [lb. /in. 3] 

The  technique  employed  in  step  b.  to  predict  the  new  t  is  described  below.  The  new  t 
is  predicted  by  passing  a  parabola  through  three  points  on  a  plot  of  t  versus  margin 
of  safety,  MS.  The  points  are  (TBAR  -0,  MS  -  -1),  (tj,  MS  tj),  (tj  +  x,  MS  (tj  + 1), 
see  Figure  4-1 7  The  new  t  is  found  by  solving  for  the  proper  root  of  the  resulting 
equation.  The  process  is  started  by  assuming  the  slope  at  t*0  to  be  0  for  the  first 
iteration. 
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Figure  4~  17  Resizing  Procedure 


.ELEMENT  OPTIMIZATION  PROCEDURE.  The  element  optimization  procedure 
is  used  to  adjust  detail  dimensions  of  an  element  so  as  to  make  the  most  efficient  use 
of  the  material  while  maintaining  a  given  weight.  For  example,  refer  to  Figure  4-13. 


Figure  4-1&.  Typical  Panel  Construction 
given  weight  within  the  boundaries  imposed 


The  panel  element  shown  contains  eight 
detail  dimensions  or  design  variables. 

All  of  the  design  variables  or  as  few  as 
one  may  be  active,  depending  on  the  choice 
of  the  user.  The  inactive  design  variables 
are  treated  as  constants  and  remain  un¬ 
changed.  The  object  of  the  optimization 
procedure  is  to  find  the  optimum  set  of 
active  design  variables,  i.e. ,  the  set  that 
describes  the  most  efficient  panel  of  a 
jy  upper  and  lower  limits  on  the  detail 


dimensions.  For  example,  suppose  that  all  of  the  design  variables  are  active  for  the 


panel  shown  in  Figure  4-18-  The  panel  weight  is  given  by: 


WP  “  (psk  *  71  *  pat  ’  Ast/B1)  *  PBBrt  Area 
where, 

-  T2  •  B2  +  T3  •  B3  +  T4  •  B4  [in. 2] 
Panel  Area  *  Panel  Length  •  Panel  Width 


Since  the  panel  area  1s  not  changing,  the  panel  weight  la  constant  If  the  unit  weight  w 
is  constant:  1 


w  *  W  /  Panel  Area 
P  P 

w  -  p  .  *  T1  +  P '  (T2  •  B2  +  T3  •  B3  +  T4  •  B4)/B1 
p  sic  St 

Given  any  seven  design  variables,  the  density  of  the  skin,  Psj.f  and  the  stiffener, 
and  a  unit  weight  wp,  the  eighth  design  variable  can  be  found  from 

T1  -  lwp  "  P8t  <T2  •  B2  +  T3  •  B3  +  T4  •  B4)/Bl]/pslc 

For  this  example,  there  are  seven  Independent  active  design  variables  and  one  depen¬ 
dent  active  design  variable,  Tl.  A  non-linear  math  programming  routine  is  employed 
to  find  the  optimum  set  of  independent  active  design  variables.  The  most  efficient 
design  is  defined  by  that  set  of  active  design  variables  for  which  the  following  func¬ 
tion  is  minimum. 

L  J 

P~£  £  F  (M^j)  .  F  (MC) 

1-1  J«1 

where:  1  denotes  the  failure  mode 

j  denotes  the  loading  condition 
MS  -  margin  of  safety 

MC  •  side  constraint  margin  for  the  dependent  active  design  variable, 
ffl~Tlmln)/Tlni)w 

The  function  F  is  given  by: 

F(x)  -  1/x  x  a  c 

F(x)  •  I x/c  (x/c  -  3)  +  3 ]/«  x  <  t 

where  c  -  0. 01. 

The  independent  active  design  variables  are  kept  within  upper  and  lower  boundaries  by 
means  of  a  mapping  function.  These  detail  dimension  variables  are  transformed  into 
optimization  variables,  and  back  into  detail  dimension  variables  by  the  mapping  function 
presented  in  Figure  4*19.  This  mapping  technique  provides  for  bounded  design 
variables  without  Imposing  constraints  on  the  optimization  variables. 

The  optimization  problem  posed  In  the  following  manner  is  an  unconstrained  ncn-linear 
mathematical  programing  problem.  A  number  of  techniques  are  available  to  solve  this 
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Figure  4-19.  Design  Variable  Mapping  Function 

problem;  the  method  used  by  APAS  in  is  the  Fletcher~Powell-Davidon  unconstrained 
minimization  technique  (Reference  7) . 

STRUCTURAL  ELEMENT  SYMMETRY  GROUPS 

A  symmetry  group  Is  a  group  of  structural  elements  that  have  Identical  designs.  When 
a  number  of  structural  elements  are  placed  into  the  same  symmetry  group,  only  one 
design  is  produced.  The  design  of  the  element  respects  all  of  the  margins  of  safety 
for  all  of  the  elements  in  the  group.  This  technique  provides  a  means  by  which 
fuselage  centerplane  symmetry  can  be  respected  without  duplicating  reversible  loading 
conditions.  It  is  often  desirable  to  make  adjacent  panel  elements  identical  for  ease  in 
manufacturing.  This  can  be  accomplished  with  symmetry  grouping  also.  Since  the 
uae  of  symmetry  groups  reduces  the  number  of  independent  design  variables  of  the 
structure,  it  can  be  significant  in  reducing  execution  time  and  should  be  employed 
wherever  possible. 

STRUCTURAL  ANALYSIS 


This  section  presents  the  techniques  used  to  calculate  the  applied  stresses,  including 
the  fatigue  stress  spectrum,  and  the  methods  used  to  calculate  margins  of  safety  for 
static  strength,  fatigue,  flaw  growth,  and  residual  strength  criteria. 

INTERNAL  LOADS  ANALYSIS.  The  internal  loads  analysis  is  based  on 
classical  box  beam  theory  (Reference  5) .  The  assumptions  made  are:  plane  sections 
remain  plane  under  the  action  of  bending  moments  and  axial  loads,  cross  sections  are 
free  to  warp  when  torque  is  applied,  and  the  structure  obeys  a  linear  elastic  stress- 
strain  law. 
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Axial  Stress.  The  axial  stresses  are  made  up  of  stresses  due  to  axial  loads, 

and  stresses  due  to  bending  moments.  Hie  equation  used  to  calculate  the  axial  stresses 

is: 

a-  EC 

where:  E  -  Modulus  of  elasticity  [psi] 

c  *  Strain  (computed  as  shown  below)  [In.  /la.  ] 

Mx  (EI^)  -  Mz  (EI^)  _  p 

C  - —  (X-X)  - J-  (z-z)  *  — - 

(EIXX)(EIU)  -  (El**)  (EIxxXEIzz)  -  (Elxj) 

where:  Mx  ■  Net  bending  moment  about  a  horizontal  axis  passing  through  the 
centroid  [in.  -lb.  ] 

Mz  ■  Net  bending  moment  about  a  vertical  axis  passing  through  the 
centroid  [in. -lb.] 

P  ■  Axial  load  [lb.  ] 

x,  z  ■  Coordinates  of  the  element,  see  Figure  4-20  [In.  ] 
x,~z  *  Coordinates  of  the  centroid,  see  Figure  4-20  fin.  ] 

Elxz.  EIxx»  Elzz*  AE  *  Section  properties,  see  •  next  section. 


Figure  4-20  Typical  Wing  Section 
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Shear  Stresses.  The  shear  stresses  resist  the  external  shear  forces  and 
the  torque  applied  to  the  section.  Under  the  basic  assumptions  the  shear  flow  is 
calculated  using  a  VQ/I  distribution  for  the  shear  forces.  The  resultant  applied 
torsion  is  due  to  the  applied  torque,  TOR,  and  the  couples  resulting  from  shifting 
the  shear  forces  XS  and  ZS  to  the  shear  center.  This  net  torsion  is  resisted  internally 
by  a  shear  flow  distributed  according  to  a  T/2A  distribution.  In  the  case  of  multiple 
cell  structures .  such  as  multi-spar  wings,  the  cells  are  assumed  to  have  equal  twist¬ 
ing  angles.  For  a  further  description  of  the  method  see  Paragraphs  17. 9  through 
17. 11  of  Reference  8. 

SECTION  PROPERTIES.  Section  properties  of  the  cross  section  of  the  wing  or 
fuselage  are  calculated  at  each  station  where  structural  sizing  is  performed.  These 
properties  are  used  to  calculate  the  Internal  loads  distribution  and  to  provide  stiffness 
information,  hi  order  to  simplify  the  calculation  of  section  properties  the  following 
assumptions  are  made:  (1)  the  material  that  resists  bending  moments  is  assumed  to 
be  smeared  uniformly  between  nodal  points,  {2)  only  tho  slcja  and  shear  webs  are 
effective  for  resisting  shear  loads  and  torsion.  The  following  equations  are  used  to 
calculate  section  properties  (see  Figure  4-20). 

EA- y*E  da 

X  •  y*£  xda 

‘  ^ 

EX,*-  Je  *2  da  -  EA •  a  * 

El,,  - y*2  da£  -  EA  •  X2 
EZxx  •  /E  xsda  -  EA  •  x  *  z 

where:  x,  z  axe  the  coordinates  of  the  incremental  area  da  [In.  ] 

EA  la  die  axial  stiffens! s  of  the  croas-aectlon  [lb.  ] 

Iif  are  the  coordinates  of  the  centroid  of  the  cross-section  [la.  ] 

EZjqc  is  the  moment  of  Inertia  of  A  taken  about  an  x  axis  passing  through  the 
centroid  multiplied  by  E  [lb. -in.2] 

EIj2  1*  moment  of  Inertia  of  A  taken  about  a  z  axis  passing  through  the 
centroid  multiplied  by  E  [lb.  -in. 2] 
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EIjjj  Is  the  product  of  Inertia  of  A  with  respect  to  the  centroid  multiplied  by  £ 

STATIC  STRENGTH  ANALYSIS.  The  depth  of  the  static  strength  analyses 
performed  by  APAS  IV  ii  consistent  with  typical  pre-design  stress  analyses.  The 
analytical  techniques  and  their  sources  are  described  in  this  section.  The  failure 
nodes  included  are  summarized  in  Table  4-18.  APAS  IV  computes  margins  of  safety 
for  each  failure  mode  and  uses  these  values  to  direct  structural  design  optimization. 
The  critical  failure  mode  margin  of  safety  is  included  in  the  computer  output  for  each 
element  and  load  condition. 


Table  4-18.  Panel  Element  Failure  Modes 


Failure  Mode 


Panel  Construction  Type 
(see  Figure  2-4) 


IDSQBQODQEn 


Local  Buckling  • 
Diagonal  Tension 


Crippling* 


Inter-Rivet  Buckling  and  Wrinkling 
Panel  General  Instability . 


Wide  Column  Buckling  *  . 

General  Yielding . 

Distortion  Energy  Theory 


Maximum  Fiber  Strain  in  a  , 


Local  Buckling  (Compression  and  Shear) 


Critical  local  buckling  stresses  for  compression  and  shear  loading  are  computed  by 
APAS  IV  using  the  following  equations.  These  equations  were  obtained  from  Reference  5 


Compression  Buckling 

F  Af 

Cr  12(l-oe2)  b 

where:  F^  -  critical  compression  buckling  stress  [psi] 

kc  -  compression  buckling  coefficient  [dimensionless] 


1 


E  -  modulus  of  elasticity  [psi] 

Ve  -  elastic  Poisson's  ratio  [dimensionless] 

t  -  thickness  [in.] 

b  -  short  dimension  of  plate  or  loaded  edge  [in.  ] 


Shear  Buckling 


F 

scr 


n2k  E 
s 

12  (1  -  ue2 ) 


where:  F 

scr 

it* 


critical  shear  buckling  stress  [psl] 
shear  buckling  coefficient  [dimensionless] 


b  -  short  dimension  of  plate  [In.  ] 

The  buckling  coefficient  is  dependent  on  the  aspect  ratio  of  panel  length/width 
and  panel  edge  fixity.  For  APAS  IV,  the  aspect  ratio  Is  assumed  to  be  large 
and  die  corresponding  asympototic  value  of  the  buckling  coefficient  for  the 
appropriate  edge  fixity  Is  used.  Typical  values  of  shear  and  compression  buckling 
coefficients  for  various  edge  fixity  conditions  are  shown  In  Figure  4-  21.  The 
actual  coefficients  used  in  APAS  IV  for  each  available  type  of  stiffened  skin  panel 
construction  are  shown  in  Figure  4-  22. 


For  some  flight  vehicle  designs,  it  may  be  a  requirement  that  buckling  of  the  skin 
panels  is  not  permitted  up  to  a  specified  percent  of  limit  load.  APAS  IV  has  the 
capability  to  handle  this  design  criterion.  Both  shear  and  compression  buckling 
are  considered.  The  interaction  equation  used  in  APAS  IV  to  combine  the  effects  of 
shear  and  compression  was  taken  from  Reference  5  and  is  presented  below: 

Rc  +  R,2  -  h  0 

where:  Rc  -  applied  compression  stress/compression  buckling  stress 
Ra  -  applied  shear  stress/shear  buckling  stress 
The  associated  margin  of  safety  equation  is: 


M.S. 


1 


l 


4-52 


k,.- 0.43 


kg  -  1.28 


kg  -  4.00 
k,-  5.40 


kg- 07 
kj-  9.00 


LEGEND:  F  -  FREE  EDGE 

SS-  SIMPLY  SUPPORTED  EDGE 
C  -  CLAMPED 


Figure  4*21.  Shear  (ks)  and  Compression  (kg)  Buckling  Coefficients  for 
Various  Edge  Fixities 


Diagonal  Tension  Analysis 

Maximum  allowable  panel  shear  stresses  are  determined  by  APAS  IV  using  the 
relationship  shown  in  Figure  4-23 .  Parameters  Fs  and  Fsu  are  the  maximum 
allowable  shear  and  the  material  ultimate  shear,  respectively.  Parameter  FSCr  is 
critical  shear  stress  at  which  shear  buckling  initiates.  The  equation  used  by  APAS 
IV  to  compute  Fscr  is  described  In  the  local  buckling  failure  mode  section. 

Crippling 


The  method  for  the  crippling  analysis  used  in  APAS  IV  was  taken  from  Reference  9. 
The  crippling  strain  for  the  combination  of  stiffener  and  effective  skin  is  computed 
by  the  following  equation: 


lcc 


r  bn  ln  fcco 

Zb  t  E 
non 


4-53 


TYPE  1  INTEGRAL  8 LAOS 


TYPE  S  EXTRUOED CLOSED  HAT 


'■A, 


Figure  4-23.  Diagonal  Tension  Chart 


where;  <cc 

En 

*ccn 


crippling  strain  for  section  [in. /in.  ] 
modulus  of  elasticity  of  element  n  [psi.  ] 
effective  element  width  [in.  ] 
element  thickness  [in.  ] 
element  crippling  stress  [psl.  ] 


The  element  crippling  stress  (fccn)  is  obtained  from  Figure  4-24. 
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THIS  CURVE  APPLICABLE  TO  ALL 
-  OUCTILE  ALLOYS  AT  BOTH  ROOM 
ANO  ELEVATED  TEMPERATURES 


Inter-Rivet  Budding  and  Sheet  Wrinkling 

Inter-rivet  buckling  involves  a  failure  of  the  skin  between  rivets.  When  both  the  skin 
and  the  stiffener  fail,  it  is  known  as  a  wrinkling  failure.  APAS  IV  checks  for  both  of 
these  failure  modes  using  the  methods  taken  from  Reference  5.  Inter- rivet  buckling 
strain  is  computed  using  the  following  equation: 


«ir“ 


£  "2  Oa]2 

i2(i-ue2)  \pJ 


where:  fjr 
C 

‘s 

P 

E 


-  inter-rivet  buckling  strain  [In. /La.] 

-  end  fixity  (C  equals  4,  for  all  cases)  [dimensionless] 

-  skin  thickness  [In.  ] 

-  rivet  pitch  [  in.  ] 

-  modulus  of  elasticity  [psl] 

-  Poisson's  ratio  [dimensionless] 


Rivet  pitch  spacing  is  set  equal  to  four  times  the  rivet  diameter  for  all  cases.  The 
rivet  diameter  for  each  case  Is  selected  based  on  skin  thickness  according  to  Table 

4-19. 


Table  4-19.  Rivet  Diameter  Versus  Skin  Thickness 
cm,  (inches) 

Skin  Thickness  (ts)  Rivet  Diameter 


0.  000  (0.  000) 

0. 064  (0. 025) 

0. 127  (0. 050) 

0.318  (0.125) 

0. 635  (0. 250) 


0.318 

(0. 1250) 

0. 397 

(0. 1563) 

0.476 

(0. 1875) 

0.  635 

(0.2500) 

Sheet  wrinkling  strain  is  computed  using  the  equation  presented  below: 

,  k.  *2  r.f 

<w  -  2*  lb”  ] 

12(l-ve)  \V 
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t 


where:  <w  -  wrinkling  strain  [in. /in.] 

kw  -  wrinkling  coefficient  [dimensionless] 
t8  -  skin  thlctaess  [  in.  1 
bs  -  stringer  spacing  [in.  ] 

The  empirical  wrinkling  coefficient  (kw)  is  a  function  of  the  effective  rivet  offset  and 
local  geometry.  The  effective  rivet  offset  is  determined  using  Figure  -*-25  and  Is 
used  in  Figure  4-26-  to  evaluate  the  wrinkling  coefficient. 


Figure  4-25.  Experimentally  Determined  Figure  4-26.  Experimentally  Determined 
Values  of  Effective  Rivet  Offset  Coefficients  for  Failure  in  Wrinkling 

Mode 


Panel  General  Instabilit 


This  section  Is  a  description  of  the  analysis  uaed  to  calculate  general  Instability 
allowables  for  panel  types  10,  11,  and  12  (see  Figure  4-8.  This  analysis  procedure 
was  taken  from  Reference  10.  Design  formulas  are  used  to  provide  conservative 
estimates  of  the  buckling  allowables. 


The  moments  of  inertia  and  stiffnesses  in  both  directions  are  calculated  for  a  plate  or 
sandwich  panel  In  the  conventional  manner.  See  Figure  4-27  for  sign  convention  used 
in  the  development  of  the  design  equations  for  buckling. 


D11  *  Ex  V(1  "  ’’xy 

°22  *  %  Vd  ~  »xy  Uyx> 
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Di2  =  Du  uyx 
DS6  “  Gxy  rx 

The  shear  buckling  of  a  simply  supported  orthotropic  plate  can  be  reasonably  estimated 
with  the  following  formula 


* I - —  2 

Nxycr  "  c3  yDll  D22  /  b 

4 _ 

where:  Ci*(b/a)  yDU/D22 

C2  ■  (D12  *  2Dg6)  /  yDji  D22 

C3  «  32. 8  +  20  C2  +  14.2  (Cx)2*4  +  24. 8  C2  Cx2 

No  simple  correction  for  the  effect  of  curvature  on  the  shear  buckling  allowable  is 
available  at  present. 

The  buckling  allowable  for  a  flat,  simply-supported ,  orthotropic  plate  under  biaxial 
loads  may  be  found  using  the  following  method.  Given  a  ratio 


or  *  Nx/Ny 

the  allowable  in  the  axial  direction  may  be  expressed  by 


^N*cr\»nel 


2  (Djjj  ♦  2D66)(n/b)2  +  D22  (n/b)4(a/m)2] 

2  2 
1  +  or  (an)  /  (bm) 


where  m  and  n  are  possible  half-wave  numbers  into  which  the  panel  may  buckle  in  the 
x  and  y  directions,  respectively.  This  formula  is  evaluated  for  the  first  five  modes 
in  each  direction,  and  the  minimum  value  is  chosen.  The  allowable  Nycr  is  obtained 
similarly. 


An  estimate  for  the  correction  due  to  curvature  on  the  compressive  buckling  allowable 
is  obtained  as  described  in  Reference  8.  The  buckling  allowable  of  the  full  cylinder 
from  which  the  panel  was  cut  is  added  to  the  flat  plate  allowable  as  obtained  above. 

For  an  orthotropic  cylinder,  the  cylinder  buckling  allowable  is  approximated  as 


(Njecr)cylinder  *  2  VExEy  1  /  CR^/3  U  '  ”xy  ^yx)  ^ 
where  t  and  R  are  the  thickness  and  radius  of  the  cylinder,  respectively. 
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The  ratios  of  the  applied  loads  to  the  allowables  are  formed: 

Rs  *  Nxy/Nxycr 

WScr 

The  buckling  margin  of  safety  for  each  panel  Is  calculated  from  the  interaction  equation 


The  above  analysis  is  brief  and  offers  various  degrees  of  approximation,  depending  on 
the  complexity  of  the  section  being  analyzed.  For  construction  with  orthotropic  flat 
panels,  it  provides  excellent  estimates  for  the  buckling  allowables  of  simply  supported 
panels.  With  the  addition  of  curvature,  it  provides  somewhat  less  accurate  allowables. 

Wide  Column  Buckling 

This  section  is  a  description  of  the  analysis  used  to  calculate  the  general  instability 
allowables  for  panel  types  1  through  9,  (see  Figure  4-8). 


Wide  column  buckling  analysis  of  multi-rib  structures  assumes  that  the  cover  panel 
behaves  as  a  simply  supported  column.  The  ribs,  oriented  perpendicular  to  the  load, 
are  assumed  to  provide  the  continuous  simple  supports.  The  effect  of  spar  support 
at  the  unloaded  edges  of  the  column  Is  ignored  In  this  analysis.  The  method  used  by 
APAS  IV  for  wide  column  analysis  was  taken  from  Reference  5  and  Is  described  below. 

The  relationship  between  critical  column  strain  versus  slenderness  ratio  (L'/p)  Is 
shown  in  Figure  4-28. 

For  large  values  of  slenderness  ratio,  a  form  of  the  Euler  column  equation  applies: 


( 

c 


2 

IT 


(L'/p  )2 


where:  e£  -  column  failing  strain  [in. An.] 


(L'/p)-  slenderness  ratio  (effective  length/ radius  of  gyration) 


r 


TRANSITION  RANGE 


Figure  4-28-  Critical  Column  Strain  Versus  Slenderness  Ratio  (L’/p) 


For  small  Values  of  slenderness  ratio,  the  critical  column  strain  transitions  from  the 
crippling  strain  to  the  Euler  critical  column  strain.  The  following  parabolic  approxi¬ 
mation  is  used  to  represent  this  transition. 


where: 


-  column  critical  strain  [In.  /In.  ] 

-  crippling  strain  [in.  /in.  ] 

-  buckling  strain  for  column  cross-section  [in.  /in.  ] 

-  Euler  column  strain  [In.  /in.  ] 


The  equation  applies  for  c  >  e  .  Yield  strain  («„)  Is  substituted  for  *  ,  when 

c  cr  y  cr 

*cr >  €y* 


General  Yielding 

To  ensure  that  elastic  stress  conditions  exist  up  to  limit  load  for  each  structural 
design,  APAS  IV  compares  element  tensile  or  compressive  stresses  to  material  yield 
•for  all  loading  conditions. 


Distortion  Energy  Theory 

The  distortion  energy  theory  (Hencky  -  Von  Mlses  theory,  Reference  11)  is  another 
failure  mode  criterion  used  by  APAS  IV.  This  theory  Is  based  on  the  assumption 
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that  failure  occurs  when  the  distortion  energy  corresponding  to  the  principal  stress 
components  equals  the  distortion  energy  at  failure  for  the  maximum  allowable  aerial 
stress.  This  failure  criterion  is  defined  by  the  equation: 


1 


al  °2  +  c22“an«: 


The  boundary  curve  defined  by  this  equation  for  all  possible  combinations  of 
principal  stresses  is  shown  in  Figure 4-29..  Any  principal  stress  combination  that 
falls  outside  this  boundary’  curve  represents  a  negative  margin  of  safety. 


Figure  4-29  Distortion  Energy  Theory 
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Laminate  Analysts 


The  composite  panels  (construction  type  12)  are  specially  orthotropic;  if  the  panel  Is 
of  sandwich  construction.  It  Is  assumed  that  the  core  supplies  no  in-plane  stiffness  but 
Is  perfectly  rigid  in  the  out  of  plane  direction,  and  that  each  of  the  composite  faces 
carries  one-half  of  the  applied  in-plane  loads. 

The  laminate  analysis  Is  designed  to  And  the  strains  in  the  0* ,  90* ,  and  *  a*  plies. 

The  laminate  longitudinal,  transverse  and  shear  strains,  <x,  fy,  and  c^y  respectively, 
are  calculated  using  the  laminate  in-plane  constitutive  matrix,  £a],  and  the  applied 
running  loads,  Nx,  Ny,  and  Nxy,  as  shown  in  Figure  4-27 . 


<x 

• 

Nx 

<y 

«  [A-1] 

Ny 

€xy 

1 

**xy 

The  laminate  strains  are  rotated  using  the  transformation  matrix.  [  T  ] ,  for  each  ply 
angle  in  the  laminate.  These  strains  are  used  for  computing  the  margins  of  safety. 

The  margins  of  safety  for  failure  of  a  laminate  of  orthotropic  materials  are  computed 
by  using  the  six  allowable  failure  strains  of  the  basic  lamina  material  and  the  orienta¬ 
tion  angle  of  each  ply  In  the  laminate.  The  strains  are: 

+  fi]_,  tension  in  the  II  direction 
+  €22*  tension  In  the  22  direction 
+  positive  shear 

-  compression  in  the  11  direction 

-  <22  •  compression  in  the  22  direction 

-  €  i2 ,  negative  shear 

The  laminate  strains  are  calculated  and  then  transformed  to  coincide  with  each  ply 
material  axis  system  as  shown  in  Figure  4-30.  The  transformed  strains  are  then 
compared  with  the  appropriate  allowable  strains  and  three  margins  of  safety  are 
obtained  for  each  ply,  for  each  loading  condition. 

The  minimum  margin  from  all  of  the  plies  then  becomes  the  final  margin  for  the 
ultimate  strain  failure  mode  of  the  laminate.  The  equation  used  for  each  margin  of 
safety  Is: 


M.S. 


f  PSAL 
C 


1 
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Figure  4-30..  Coordinate  Transformation  for  the  or  Ply 

where  cPSAL  Is  the  ultimate  strain  allowable  and  c  is  the  applied  strain. 

FATIGUE  STRESS  SPECTRUM.  The  fatigue  stress  spectrum  is  based  on  the 
flight  profile  and  load  spectrum  discussed  in  previous  sections.  It  is  made  up  of  a  group  of 
minimum  and  maximum  stresses  and  the  number  of  applications  expected  during  the 
design  life.  This  spectrum  is  used  for  the  fatigue  analysis  and  the  flaw  growth 
analysis  of  the  element  discussed  in  this  section. 

Minimum  and  maximum  stresses  are  calculated  for  each  subsegment  of  the  fatigue 
spectrum  (see  Table  *1-13)  •  These  stresses  are  calculated  from  the  segment  constant 
stress  oc  and  the  segment  alternating  stress  ez. 

°min  m  ac  -  Aff*  oa 

°max  3  sc  +  A*  '  ffa 

The  value  of  -c  and  oa  **~e  in  general  different  for  each  segment  and  are  calculated  by 
forming  linear  combinations  of  the  stress  due  to  the  fatigue  spectrum  conditions  (see 
Table  4-ltf . 


where  i  denotes  the  fatigue  spectrum  condition  number  and  j  denotes  the  segment  number. 
The  constants  c..  and  a,. are  based  on  the  flight  profile  (see  Table  4-12)  and  are  stored 
within  the  program. 

he  ground-air- ground  (G-A-G)  cycle  shown  in  Figure  4-31  was  not  previously  defined. 
It  dominates  fatigue  damage  and  flaw  growth  in  many  areas  of  the  structure  of 
transport  aircraft.  This  stress  excursion  Is  due  in  part  to  the  difference  between  the 
groundborne  load  distribution  and  the  airborne  distribution,  and  in  part  to  cabin  pres¬ 
surization.  A  G-A-G  spectrum  Is  calculated  automatically  within  the  program  at  each 
analysis  point.  The  G-A-G  cycle  is  defined  as  the  maximum  stress  excursion  between 
the  peak  inflight  stress  (i.e. ,  the  maximum  gust  occurring  in  that  flight)  and  the  peak/ 
valley  groundborne  stress  (i.  e. ,  the  maximum  taxi  Ag).  Several  high  peaks  of  cyclic 
loads,  such  as  those  due  to  gust  encounters  in  stormy  weather,  tend  to  occur  on  the 
same  flight.  It  would  therefore  be  conservative  to  use  ail  peak  loads  expected  in  the 
total  aircraft  life  in  building  the  G-A-G  spectrum.  To  avoid  tMs  overconservatism  a 
frequency  factor  is  introduced,  which  has  the  effect  of  skipping  ever  some  of  the  peak 
loads.  A  frequency  factor  equal  to  two  is  considered  appropriate  for  transport  aircraft. 
Thus,  every  other  peak  is  included  in  the  G-A-G  spectrum.  Frequency  factor  is  a 
user  input. 


A  unique  stress  spectrum  is  generated  for  each  structural  element  based  on  the  local 
stress  history  and  is  used  for  the  fatigue  analysis  and  flaw  growth  analysis  also  reported 


Figure  4-31.  Simplified  Flight  Profile 
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in  this  section.  No  provision  is  currently  available  for  changing  the  load  profile; 
however,  additional  load  spectra  can  be  incorporated  into  the  existing  program  with 
very  little  programming  effort. 


FATIGUE  ANALYSIS.  Fatigue  damage  is  defined  as  the  ratio  of  the  number  of 
applied  stress  cycles,  n,  of  a  given  stress  magnitude  to  the  number  of  allowable  stress 
cycles,  N,  of  the  same  stress  magnitude.  Miner's  Rule  (Reference  12)  is  the  basis  of 
fatigue  damage  analysis  performed  by  the  subroutine,  PRODAM.  Under  this  concept, 
fatigue  damage  is  assumed  to  be  linearly  cumulative,  and  fatigue  failure  Is  assumed 
to  occur  when  the  damage  summation  equals  unity. 


Fatigue  Damage 


Fatigue  Failure 


n 

m 


N_ 


m 


To  facilitate  the  analysis.  S-N  curves  are  plotted  from  test  data  for  several  values  of 
stress  ratio,  R.  Allowable  cycles  for  each  subsegment  are  read  from  the  curves 
as  shown  in  Figure  4-32. 

where,  ‘  s_. 

R  »  (Also  see  Figure  2-27) 
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Figure  4-32.  Fatipie  Damage  Determination 


A  review  of  previous  General  Dynamics  programs  and  other  sources  did  not  produce 
nearly  enough  component  S-N  data  to  fill  the  required  data  bank  indicated  by  Table  4 -21). 
Much  of  the  data  reviewed  was  generated  for  specific  configurations  and  load  spectra. 
Manufacturers  usually  test  splices  and  other  fatigue  critical  details  but  seldom  develop 
S-N  curves  for  typical  structure  and  spectra.  Even  less  data  is  published  because  com¬ 
ponent  test  results  are  frequently  considered  proprietary  or  sensitive  to  a  particular 
project. 


Table  4  -20.  Availability  of  Fatigue  Data 
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In  some  cases  the  component  data  was  incomplete.  To  facilitate  extrapolation,  curves 
of  stress  versus  stress  ratic  at  constant  cycle  values  were  plotted  from  the  original  data. 
Expanded  S-N  curves  were  . — m  drawn  based  on  the  extrapolated  data.  For  the  many 
cases  where  component  data  was  not  available,  reduction  factors  were  applied  to  un¬ 
notched  coupon  data  for  the  appropriate  material.  A  complete  set  of  data  was  gener¬ 
ated  by  this  method.  However,  S-N  curves  plotted  from  this  data  did  not  show  the 
trends  and  consistency  expected.  Anomalies  in  the  component  data  due  to  inconsistent 
test  parameters  were  still  present  in  the  expanded  S-N  curves.  A  complete  and  con¬ 
sistent  set  of  S-N  curves  for  all  required  component  types  could  not  be  obtained  with 
this  approach. 

Subsequently,  a  second  method,  Reference  11,  for  plotting  S-N  curves  from  limited 
data  was  employed.  This  procedure  utilizes  two  Hewlett-Packard  9100B  computer  pro¬ 
grams  to  curve  fit  and  plot  data.  More  reliance  is  placed  on  un-no tched  coupon  aata, 
and  component  fatigue  notch  factors  are  plotted  versus  life  to  ensure  that  smooth  and 
consistent  S-N  curves  are  generated. 
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Data  from  constant  life  cuts  of  these  curves  is  stored  in  the  program.  An  interpolation 
routine  is  used  in  the  program  to  retrieve  allowable  cycles  from  the  stored  data. 

This  rather  simple  approach  is  widely  used  in  fatigue  life  predictions  of  transport 
aircraft.  The  more  severe  load  spectra  of  fighter  type  aircraft  produce  more 
significant  residual  stresses  at  points  of  stress  concentration  and  may  warrant  a 
more  sophisticated  analytical  treatment. 

CRACK  GROWTH  ANALYSIS 

The  crack  growth  analysis  procedure  used  in  the  APAS  IV  Program  predicts  how 
a  crack  grows  under  cyclic  fatigue  spectrum  loading.  Two  basic  rate  equations 
are  used,  one  for  load  steps  with  positive  stress  ratios  and  the  other  for  negative 
stress  ratio  load  steps.  The  Walker  effective  stress  concept  introduced  in  the 
Paris  rate  equation  (Reference  13) ,  serves  as  the  basis  for  predicting  crack 
growth  rates  with  positive  stress  ratio  load  steps.  The  Chang  crack  growth  rate 
model  (Reference  2)  is  used  for  negative  stress  ratio  load  steps.  Load  interac¬ 
tion  effects  are  accommodated  by  employing  the  Willenborg  effective  stress  con¬ 
cept  with  the  Chang  acceleration  scheme.  The  model  is  identified  as  Willenborg/ 
Chang  model  in  Reference  2.  These  crack  growth  prediction  methods  are 
implemented  in  subroutines  FREGRO,  CORREC ,  CRITIC,  CLAMDA  and  FLTGRO. 

Paris  formulated  a  crack  growth  rate  model  as  a  function  of  the  stress  intensity 
factor  range,  AK,  and  empirical  material  dependent  constants  (Reference  14). 


da 

dN 


F  (AK,  empirical  constants) 


where  da/dN  is  the  cyclic  growth  rate  and  the  stress  intensity  factor  range,  aK,  is  a 
measure  of  the  crack  tip  stress  field  (see  Figure  4-33) 

-  K  =  A  a  0  (a)  v/rTa 

=  K  -  K  . 
max  min 

A<r  =  (r  -  a  .  ,  range  of  remotely  applied  cyclic  stresses  [psil 
max  min 

a,  half  crack  length  [in.  ] 

0(a) ,  correction  factor  that  accounts  for  geometric  effects  [dimensionless] 


Subsequent  to  the  Paris  formulations,  results  of  constant  amplitude  tests  identified 
variations  in  crack  growth  rates  due  to  different  stress  ratios.  Based  only  on  positive 
stress  ratios,  Walker  introduced  the  concept  of  an  effective  stress,  <r . 

f  <1-k>“'',nax  Rs0 

and 

R  —  / 

min  max 

m  =  stress  ratio  collapsing  constant 
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Figure  4-33.  Fatigue  Crack  Loading 


The  following  rate  equation  is  obtained  by  using  the  Paris  model  with  the  Walker 
effective  stress. 

=  F(AK,  empirical  constants) 

=  C(AK)n 

=  C[(l  -R)mamax/3(a)^a]n 

=  C[(l  -  R)m  K  ]n 
maxJ 

=  C[(l  -  R)m_1AK]n  R  >  0 

where  C  and  n  are  empirical  constants  determined  from  constant  amplitude  tests. 

Based  on  the  observation  of  the  differences  between  constant  amplitude  tension- 
compression  cyclic  test  results  and  R  =  0  data,  in  Reference  2,  Chang  proposed 
a  rate  equation  which  correlates  crack  growth  behavior  due  to  negative  stress 
ratio  cyclic  loading.  The  equation  form  is: 

dN  =  Ci(l.R2)qKmaJt]n  ,  R<0 

where  C  and  n  are  the  same  constants  as  that  for  the  positive  stress  ratio  rate 
equation  and  q  is  the  acceleration  index  which  is  determined  by  tension-compression 
cyclic  test  for  a  specific  value  of  negative  stress  ratio. 
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For  low  stress  intensity  factor  ranges,  crack  growth  rates  are  reduced  and  at  a 
threshold  stress  intensity  factor  range,  AKxH>  no  discernible  growth  is  observed. 
This  characteristic  is  accounted  for  by  considering  zero  growth  when  the  stress 
intensity  factor  range  is  equal  to  or  less  than  the  threshold  value. 


^  =  0,  AK<AKTH 


where  AKxh  =  AKxHo  (  [  -  A  [R  J )  and  AKxHo  *s  the  threshold  value  of  AK  at 
R  =  0,  A  is  a  material  dependent  constant. 

The  foregoing  rate  equations  characterize  subcritical  crack  growth  behavior  under 
constant  amplitude  cyclic  loading.  However,  application  of  these  equations  to 
variable  amplitude  spectrum  loadings  have  been  shown  to  result  in  significant 
differences  between  predicted  lives  and  actual  test  results.  The  occurrence  of 
tensile  overloads  have  been  shown  to  retard  the  crack  growth  on  subsequent  load 
steps.  Compressive  load  immediately  following  the  tensile  overload  reduces  the 
retardation  effect.  The  Willenborg /Chang  scheme  was  selected  to  account  for 
these  load  interaction  effects.  The  method  is  based  on  defining  effective  stress, 
°eff»  effective  stress  intensity  factor,  Keff,  and  the  effective  stress  ratio,  Reff. 
to  model  the  interaction  effects.  The  following  paragraphs  describe  the 
Willenborg/Chang  load  interaction  model  briefly: 


The  plane  stress  plastic  zone,  Zql>  at  the  crack  tip  due  to  a  tensile  overload  is 
given  by: 


Z 


OL 


JL 

2ir 


R, 


OL 


>0 


where  fyy  is  the  material  tensile  yield  stress  [psi] . 

When  the  tensile  overload  is  followed  by  a  compressive  load,  the  Chang  scheme 
considers  a  reduced  effective  plastic  zone  as  shown  below. 


The  plastic  zone  radius,  roL>  due  to  the  overload  step  is  depicted  in  figure  4-34. 
rOL  =  aOL  +  ZOL 


where  aoL  is  the  half  crack  length  following  the  overload  step. 
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On  any  subsequent  load  step,  i,  crack  growth  is  considered  to  be  affected  by  the 
residual  plastic  zone  radius,  roL*  if  the  plastic  zone  for  that  step,  n,  lies  within 
that  resulting  from  the  overload  condition.  The  Willenborg  model  identified  an 
effective  stress  parameter  based  on  the  stress  required,  «req  to  break  through  the 
overload  plastic  zone. 
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^req  rOL  1  2rr  1  L 


ty 
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req 
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Gallagher  modified  the  Willenborg  retardation  model  by  introducing  a  pro¬ 
portionality  constant ,  <t> .  The  effective  stress  is  then  defined  as  a  measure  of 
the  difference  between  the  stress  required  to  penetrate  the  residual  overload 
plastic  zone  and  the  actual  spectrum  stress  (Reference  15). 


and  the  proportionality  constant,  <j> ,  is  defined  as: 


<b  =• 


rhnax 


R 


shut 


where  Rghut  shutoff  ratio,  material  dependent  constant. 

The  effective  maximum  and  minimum  stress  intensity  factors  and  corresponding 
effective  stress  ratio  are  defined  as 


(Km“)efr('Wslff  m  ^ 

M-’M-  W  ^ 
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and 


Reff  |<rmin)e^y/|(rmaxje£f  Rcut  <Reff  ~Rcut 

AK  =  K  -K.=/K  \  -  /  K  .  \ 

max  mm  I  maxjeg  I  mm J  ^ 

where 

R~ut  =  negative  stress  ratio  cutoff  value 
R+  =  positive  stress  ratio  cutoff  value 


Thus  due  to  load  interaction  effects,  the  crack  growth  rate  equations  ar< : 


The  damage  accumulation  scheme  used  to  predict  crack  growth  life  consists  of  a 
three-step  procedure  performed  by  PREGRO  and  subroutines  CLAMDA  and  FLTGRO. 
In  the  first  step,  PREGRO  uses  the  Vroman  linear  approximation  method  for  a  unit- 
block  flight  spectrum  to  obtain  crack  growth  rate  per  flight  (da/dF)j  and  a  measure 
of  the  stress  intensity  factor  Kj  for  j  values  of  initial  crack  sizs.  The  second  step 
consists  of  using  a  least-square-fit  procedure  for  the  (da/dF)j  versus  Kj  values  to 
characterize  an  equivalent  growth  per  flight  rate  equation  in  subroutine  FLTGRO  uses 
the  crack  growth  per  flight  rate  equation  to  calculate  crack  growth  life. 

The  linear-approximation  method  assumes  that  the  growth  rate  is  constant  through¬ 
out  a  load  step  in  a  spectrum  so  that  the  crack  size  is  in  a  linear  relationship  with  the 
number  of  load  cycles.  The  damage  accumulation  scheme  proceeds  by  considering  a 
load  step  (i)  and  using  crmaxi  and  crmini  to  calculate  da/dN.  The  relatively  small 
incremental  change  in  crack  length,  6  a,  of  0.  Ola  is  used  to  provide  reasonable  com¬ 
putational  accuracy  in  the  following  procedure  (Reference  16) . 
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The  value  of  (6a)/(da/dN)  is  then  compared  to  the  cycles  in  that  load  step,  Ni, 
where  "a"  is  the  crack  size.  If  (6a)/(da/dN)  is  greater  than  Ni,  then  the  crack  growth 
for  that  particular  load  step  is  Aa  =  Ni  x  (da/dN);  "a"  is  increased  by  Aa,  and  the  pro¬ 
gram  proceeds  to  the  next  load  step. 

If  (6a)/(da/dN)  is  less  than  or  equal  to  Ni,  then  the  number  of  cycles  to  grow  (6a) 
is  (6a)/ (da/dN).  This  value  is  subtracted  from  Ni,  the  crack  size  "a"  is  increased  by 
(  a),  and  the  load  step  is  reconsidered.  This  process  continues  with  (6a)/(da/dN) 
being  compared  to  the  remaining  cycles  in  the  step.  The  process  is  repeated  until  all 
load  steps  in  the  block  (or  flight)  are  exhausted. 

The  cycle-by-cycle  crack  growth  analysis  is  performed  for  a  unitblock  load  spectrum 
for  j  values  of  initial  crack  size.  The  crack  growth  rate  resulting  from  the  unitblock 
spectrum  defined  for  a  period  of  Na  flights  is  then  equal  to  the  crack  growth  Aaj 
divided  by  Na. 


Aa. 


The  second  step  in  the  computational  procedure  consists  of  characterizing  the  flaw 
growth  for  the  complex  flight  spectra  into  an  equivalent  constant- amplitude  loading 
that  will  produce  the  same  crack  growth  life.  This  method  is  based  on  the  observa¬ 
tion  of  crack  growth  rate  per  flight  as  a  function  of  a  measure  of  the  stress  intensity 
factor,  K,  representing  the  unitblock  spectrum.  This  relationship  is: 


§  -  C(K)X 

and 

k  =p^)1/2  m 

where  (A<x2)*^  js  the  root -mean -square  of  the  stress  range  history. 

The  power  exponent  \  and  the  growth  rate  constant  C  are  calculated  in  sub¬ 
routine  CLAMDA  by  applying  a  least -square -fit  procedure  to  the  log  (da/dF)j 
versus  log  (Kj)  data  plot. 

In  the  third  step,  subroutine  FLTGRO  uses  the  linear-appr<  dmation  method 
on  the  crack  growth  per  flight  rate  equation  to  determine  the  crack  growth  over 
the  prescribed  crack  length  interval. 
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The  foregoing  procedure  is  applicable  to  any  flaw  geometry  for  which  the  stress 
intensity  correction  factors,  /3,  are  known.  The  program  currently  contains  factors 
for  a  wide  range  of  stiffened  panels  with  through  cracks.  Curves  for  L(a)  and  /3(a)  for 
the  case  of  a  crack  extending  equally  on  both  sides  of  a  riveted  stiffener  (illustrated 
in  Figure  4-35  are  stored  within  the  program  in  the  form  of  data  tables.  For  further 
information  concerning  the  derivation  of  these  curves  the  reader  is  referred  to 
Reference  17. 


Figure  4-35.  Stiffened  panel  crack  geometry. 

The  program  currently  contains  75  sets  of  data  for  L(a)  and  /3(a)  covering  a  wide 
range  of  stiffener  spacing  and  percent  stiffening,  including  cases  for  broken  stiffeners. 
Figure  4-36  presents  a  typical  set  of  curves.  Linear  interpolation  is  used  to  deter¬ 
mine  L(a)  and  /3(a)  curves  for  cases  that  lie  between  data  sets.  These  curves  are  used 
for  all  riveted-stiffener  plate  combinations,  (e.g. ,  panel  types  4  through  9  of 
Figure  4-4). 

For  the  case  of  integral  construction  (e.g. ,  panel  types  1,  2,  and  3),  the  panel  is 
treated  as  a  flat  plate  without  stiffeners  with  a  thickness  equal  to  t  (i.  e. ,  /3(a)  =  1. 0). 
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Figure  4-36.  Stiffened  panel  stress  intensity  correction  factors. 

For  all  stiffened  panel  types  (types  1  through  9),  panel  half  width  is  assumed  to  be 
equal  to  10  stiffener  spacings,  and  the  stress  intensity  correction  is  revised  to  include 
the  width  correction  factor. 


/3(a)  =  L(a) 


where  w  is  the  panel  half  width  [in.  ] 


On  plate  construction  types  (types  10  and  11),  no  adjustment  is  made  for  panel  width 
and  0(a)  =  1.  0.  For  plate  concepts  crack  life  is  based  on  crack  growth  from  an  initial 
length  to  a  critical  flaw  size. 


RESIDUAL  STRENGTH  ANALYSIS 

The  residual  strength  analysis  determines  the  failing  strength  of  a  damaged  panel 
and  is  performed  by  subroutine  RESID.  Damage  consists  of  skin  cracks  and  broken 
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stiffeners.  The  residual  strength  of  a  damaged  panel  is  defined  as  the  maximum  stress 
level  that  can  be  applied  to  the  panel  without  the  crack  growing  unstably  to  failure. 
Unstable  crack  growth  occurs  when  the  applied  stress  intensity  factor,  K,  exceeds  the 
fracture  toughness  of  the  skin  material,  Kc  ■ 

Unstable  crack  growth  is  allowed  to  occur  at  stress  levels  below  the  residual  strength 
of  a  panel  as  long  as  the  crack  growth  eventually  arrests  at  a  larger  crack  size.  When¬ 
ever  stress  level  of  the  most  highly  loaded  stiffener  exceeds  the  ultimate  tensile  strength 
of  the  stiffener,  it  fails,  and  the  applied  stress  intensity  factors  of  the  skin  are  recal¬ 
culated  to  reflect  the  broken  stiffener. 

Figure  4-37  illustrates  a  typical  example  of  the  residual  strength  analysis  procedure. 
The  curves  shown  are  generated  by  calculating  the  gross  panel  stress  that  causes 
stiffener  failure  and  the  gross  panel  stress  that  cause  unstable  crack  growth.  The 

following  equations  are  used  for  these  calculations. 

« 

(sheet) 

cr(unstable  crack  growth)  v,Tra 

F^u  (stiffener) 

acg  (stiffener  failure)  L(a) 


UUTUL  BMUSt »  MOM*.  •  •  •, 

nan 


Figure  4-37.  Typical  example  of  residual  strength  analysis. 
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4.1.4  DECK  STRUCTURE  .  The  input  data  deck  for  APAS  is  described  in  this 
section.  The  overall  deck  setup  is  shown  in  Figure  4-38. 
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INPUT) 


Figure  4-38.  APAS  Input  Deck  Set-up. 


4.1.5  RELATIONSHIP  OF  INPUT  TO  OUTPUT.  A  representative  input  data 
listing  for  the  APAS  program  is  shown  in  Figure  4-39.  The  input  is  entered  through 
the  NAMELIST  capability  as  discussed  in  sections  4. 1. 1  and  4.1.2.  The  results 
of  the  analysis  are  printed  out  (see  Section  4.  2)  and  also  stored  on  a  series  of  output 
files  which  may  be  later  accessed  by  parts  prediction  and  cost  routines. 
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Figure  4-J9.  A  FAS  Input  Listing 
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4.2 


RESULTS  OF  OPERATIONS 


4.2.1  DESCRIPTION  OF  RESULTS.  The  results  of  the  APAS  analysis  are  provided 
as  printed  output  and  also  stored  on  files  to  be  available  for  later  use  by  parts 
prediction  and  cost  routines. 

4.2.2  OUTPUT  FORMAT  AND  CONTENT.  A  typical  print-out  is  illustrated  in 
Figure  4-46.  The  test  case  evaluates  an  all-metal  wing  with  2024-  762  aluminum  skins, 
riveted  7075- T6  aluminum  "J"  stiffeners,  and  6A1-4V  titanium  ribs  and  spare. 

The  example  includes  fatigue,  flow  growth,  and  residual  strength  analysis. 
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SECTION  V 

OPERATING  INSTRUCTIONS 


5.1  INITIATE  THE  COMPUTER  PROGRAM 

This  computer  program  runs  on  the  Network  Operating  System  (NOS  BE  1. 0).  It  is 
written  in  FORTRAN  IV  for  use  with  the  FTN  compiler. 

This  section  shows  a  set  of  control  cards,  Figure  5-1,  (catalog  procedures)  for  oper¬ 
ating  this  program.  The  first  card  is  a  GEMS  control  card  which  puts  the  procedure 
into  the  system.  This  is  covered  in  more  detail  in  the  users  manual  volume  for  GEMS. 
This  procedure  attaches  a  preprocessor,  executes  it  to  retrieve  data  from  the  data 
base  and  puts  it  into  the  format  needed  by  the  program.  Before  attaching  and  executing 
the  program,  sense  switch  1  is  set  to  signal  the  program  that  the  data  is  coming  from 
Tape  2  instead  of  Tape  5  input.  After  the  program  has  executed,  the  post-processor 
is  attached,  is  executed  to  store  data  into  the  data  base  for  use  by  other  program 
modules.  Finally,  the  GEMS  procedure  is  called  to  bring  the  executive  back  into  core. 

This  program  can  operate  using  Tape  5  as  the  input  file  by  removing  the  switch  1  and 
pre/post  processor  control  cards.  This  then  will  operate  as  a  standalone  program, 
independent  of  the  overall  system. 

5.2  UPDATING  PROCEDURES 

This  program  can  be  mndifiari  by  using  the  standard  CDC  UPDATE  utility, 
that  is  by  using  the  *1  and  *0  for  inserting  and  deleting  desired  code.  However,  if 
the  program  is  modified,  new  relocatable  binary  and  task  files  will  need  to  be  made. 
The  catalog  procedure  shown  in  Figure  5- 1  will  accomplish  this  by  taking  the  follow¬ 
ing  steps: 

ATTACH ,  PRO  FIL,  FILENAME ,  ID  -NAME 
BEGIN, PLBDECK,  ,PL,LG,TSK,IDN. 

NOTES: 

.  Procedure  is  attached  as  PROFIL. 

.  PLBDECK  is  name  of  procedure. 

.  PL  ■  file  name  of  the  OLDPL. 

.  LG  -name  desired  for  relocatable  binary. 

.  TSK  ■  name  of  desired  TASK  FIL. 

.  IDN  -  desired  ID  name. 

For  information  involving  UPDATE  DECK  names,  refer  to  the  compilation  listings  of 
the  program.  For  Information  involving  the  use  of  the  UPDATE  capabilities,  refer  to 
the  Systems  UPDATE  Manual. 

5-1 


CATLG,PROCS.ST£F.aPfiS.  BATCH 
ATTACH, PRE.PRc POST  ,  ID= R£ SO  760 L7 . 
PRE,*GATA. 

REWIND,  APASI. 

RETUR*  ,  PRS, 

IF,SSW,-1,1» 

SWITCH, i. 

IF,SSW, -6,1. 

SWITCH, 6. 

IF, BATCH, 1. 

SWlTCh, 6. 

ATTACH,  TASK,AF  ASTSK,  IO  =  REED?6C<.7. 
TASK, 

RETURN, TASK, APASI. 

ATTACH, POST,PREPOST, 10 =®EEQ?€ 2*7. 

REWIND, APA SOI,  APAS02. 

POST »*0ATA. 

RETURN, POST, APASOi, APA SO  2. 

RETURN, *OATA. 

BEGIN, GEMS,  ACJCL, 

/OATA 

APASI 

/ECR 

APASOI 

LINKSPEC 

APAS02 

LINKSPEC 

/ECR 

/EC  F 

ENC 


Figure  5-1.  Procedure  for  Executing  APAS. 


This  program  uses  segmentation  instead  of  overlay.  Figure  5-2  represents  the 
segmentation  directives  used  for  the  segload  sequences. 
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Figure  5-2.  Segmentation  Directives 
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PLBDECK  (PL,  LG,  TSK,  ISG,  IDN) 
ATTACH  (OLDPL,  PL,  ED=IDN) 
UPDATE  (F ,  I=DUM) 

RETURN  PUM) 

RETURN  (OLDPL) 

REQUEST  (LGO,  *PF) 

FTN  (I=COMPILE,R=3) 

CATALOG  (LGO ,  LG ,  ID=IDN) 
REWIND  (LGO) 

MAP  (ON) 

ATTACH  (SEG ,  ISG ,  ID=IDN) 
SEGLOAD  (I-SEG.B-ABS) 

LOAD  (LGO) 

NOGO. 

RETURN  (LGO) 

REQUEST  (TASKFIL,  *PF) 
EDITUB  (USER,  I=PROFIL) 
CATALOG  (TASKFIL, TSK,  ID*IDN) 
RETURN  (A BS, TASKFIL) 
REVERT. 


LIBRARY  (TASKFIL, NEW) 
REWIND  (ABS) 

ADD  (*, ABS) 

FINISH. 

ENDRUN. 


Figure  5-3.  Procedure  for  Creating  a  TASKFIL 


5-4 


m 


REFERENCES 


1.  Oman,  B.  H.,  Kruse,  G.  S.,  and  Reed,  T.  F.,  "Structural  Technology 
Evaluation  Program  (STEP),  Volume  III  -  User's  Manual  for  Structural 
Synthesis,"  AFFDL-TR-77-110,  Volume  III,  Air  Force  Flight  Dynamics 
Laboratory,  WPAFB,  Ohio,  1978. 

2.  Chang,  J.  B.,  R.  M.  Hiyama,  and  Szamossi,  M.,  "Improved  Method  for 
Predicting  Spectrum  Loading  Effects  -  Final  Report,  Volume  I  -  Technical 
Summary"  AFWAL-TR-81-XXXX ,  Volume  I,  Air  Force  Wright  Aeronautical 
Laboratory,  WPAFB,  Ohio,  1981. 

3.  Tanner,  C.  J.,  Kruse,  G.  S.,  and  Oman,  B.  H.  "Computer  Program  to 
Assess  Impact  of  Fatigue  and  Fracture  Criteria  on  Weight  and  Cost  of 
Transport  Aircraft,"  (NAS1-12506) ,  NASA-CR-132648,  June  1975. 

4.  Shanley,  F.  R.,  "Weight-Strength  Analysis  of  Aircraft  Structures,"  Dover 
Publications,  Inc.,  New  York,  N.  Y.,  March  1960. 

5.  Bruhn,  E.  F.,  "Analysis  and  Design  of  Flight  Vehicle  Structures,"  Tri- 
State  Offset  Co.,  Cincinnati,  Ohio,  1965. 

6.  Anon,  Military  Specification,  "Airplane  Strength  and  Rigidity  Reliability 
Requirements,  Repeated  Loads,  and  Fatigue,"  MIL-A-8866B,  June  1975. 

7.  Fletcher,  R.,  and  Powell,  M.  J.  D.,  "A  Rapidly  Convergent  Descent 
Method  for  Minimization,"  The  Computer  Journal,  Vol.  6,  April  1963  - 
January  1964. 

8.  Peery,  D.  J.,  "Aircraft  Structures,"  McGraw-Hill  Book  Company,  New 
York,  N.Y.,  1950. 

9.  Anon,  "Astronautics  Structures  Manual,"  NASA  Marshall  Space  Flight 
Center. 

10.  Blackmon,  C.  M.,  and  Eisenmann,  J.  R.,  "Advanced  Composite  Fuselage 
Section  Optimization,"  General  Dynamics  Report  FZM-5686,  June  1976. 

11.  Eshback,  O.  W.,  "Handbook  of  Engineering  Fundamentals,"  John  Wiley  & 
Sons,  Inc.,  New  York,  May  1966. 

12.  Miner,  M.  A. ,  "Cumulative  Damage  in  Fatigue,"  Journal  of  Applied 
Mechanics ,  Volume  12,  No.  3,  Sept.  1945. 

13.  Walker,  K. ,  "The  effect  of  Stress  Ratio  During  Crack  Propagation  and 
Fatigue  for  2024-T3  and  7075-T6  Aluminum,"  ASTM  STP  462,  American 
Society  for  Testing  and  Materials,  1970. 


6-1 


Paris,  P.D.,  Gomez,  M.  P.,  and  Anderson,  W.  E.,  "A  Rational  Analytical 
Theory  of  Fatigue,"  The  Trend  in  Engineering,  Vol.  13,  No.  1,  1961. 

Gallagher,  J.  P.,  "A  Generalized  Development  of  Yield- Zone  Models," 
AFFDL-TM-74-28,  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio,  1974. 

Szamossi,  M.,  "Crack  Propagation  Analysis  by  V roman's  Model,  Program 
EFFGRO,"  NA-72-94,  Rockwell  International,  Los  Angeles,  1972. 

Poe,  C.  C.,  Jr.,  "Stress-Intensity  Factor  for  a  Cracked-Sheet  with 
Riveted  and  Uniformly  Spaced  Stringers,"  NASA-TR-358,  Washington, 
D.C.,  May  1971. 


